The ground state construction of the two-dimensional 
Hubbard model on the honeycomb lattice 
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In these lectures I consider the half-filled two-dimensional (2D) Hubbard model 
on the honeycomb lattice and I review the rigorous construction of its ground state 
properties by making use of constructive fermionic Renormalization Group methods. 

1. INTRODUCTION 

There are very few quantum interacting systems whose ground state prop- 
erties (thermodynamic functions, reduced density matrices, etc.) can be com- 
puted without appro ximations. Among these, the Luttinger and the Thirring 
model 



37, kO 



45 



46| (one-dimensional spinless relativistic fermions) , the one- 



dimensional Hubbard mode 



the Lieb-Liniger model 34 



16 



35 



36| (non-relativistic lattice fermions with spin). 



one- dimensional bosons with repulsive delta in- 



teractions) and the BCS model ll, 15, 42 1 (d-dimensional spinning fermions with 
mean field interactions, o? > 1); the construction of the ground states of these sys- 
tems is based on some remarkable exact solutions, which make use of bosonization 
techniques and Bethe ansatz. Unfortunately, in most cases these exact solutions 
do not allow to compute the long-distance decay of the n-particles correlation 
functions in a closed form (with some notable exceptions, namely the Luttinger 
and Thirring models, see 



32 



38l|). Moreover, even the computation of the thermo- 
dynamic functions crucially relies on a very special choice of the particle-particle 
interaction; as soon as the interparticle potential is slightly modified the exact 
solvability of these models is destroyed and no conclusion on the new "perturbed" 
system can be drawn from their solution. This is, of course, very annoying and 
unsatisfactory from a physical point of view. 

There is another rigorous powerful method that allows in a few cases to fully 
construct the ground state of a system of interacting particles, known as Renor- 
malization Group (RG). This method, whenever it works, has the advantage that 
it provides full information on the zero or low temperatures state of the system, 
including correlations and that, typically, it is robust under small modifications 
of the interparticle potential; even better, it usually gives a very precise mean- 
ing to what "small modifications" means: it allows one to classify perturbations 
into "relevant" and "irrelevant" and to show that the addition of a small ir- 
relevant perturbation does not change the asymptotic behavior of correlations. 
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Unfortunately, it only works in the weak coupling regime and it has only been 
successfully applied to a limited number of interacting quantum systems. Most 
of the available results on the ground (or thermal) states of interacting quan- 
tum systems obtained via RG concern one-dimensional (ID) weakly interacting 



fermions, e. g., ultraviolet 0{N) models with N > 2 |l7l . |22|| . ultraviolet QED- 



ilj: 

B 



like models 33|] and non-relativistic lattice systems 

011 

|23| . In more than one 

dimension, most of the result derived by rigorous RG techniques concern the fi- 
nite temperature properties of two-dimensional (2D) fermionic systems above the 



BCS critical temperature 



aa 



Two remarkable exceptions are 



[131, HQ. 

the Fermi liquid construction by Feldman, Knorrer and Trubowitz 191], which 
concerns zero temperature properties of a system of 2D interacting fermions with 
highly asymmetric Fermi surface, and the ground state construction of the short 
range half-filled 2D Hubbard model on the honeycomb lattice by Giuliani and 
Mastropietro 2^, l25|], which will be reviewed here. 

The 2D Hubbard model on the honeycomb lattice is a basic model for describ- 
ing graphene, a newly discovered material consisting of a one-atom thick layer 
of graphite [41|, see 12] for a comprehensive and up-to-date review of its low 
temperature properties. One of the most remarkable features of graphene is that 
at half-filling the Fermi surface is highly degenerate and it consists of just two 
isolated points. This makes the infrared properties of the system very peculiar: 
in the absence of interactions, it behaves in the same way as a system of non- 
interacting (2 + l)-dimensional Dirac particles 47|. Therefore, the interacting 
system is a sort of (2 -|- l)-dimensional QED, with some peculiar differences that 
make its study new and non-trivial 29|, |44| . 

The goal of these lectures is to give a self-contained proof of the analyticity of 
the ground state energy of the Hubbard model on the 2D honeycomb lattice at 
half filling and weak coupling via constructive RG methods. A simple extension 
of the proof of convergence of the series for the specific ground state energy 
presented below allows one to construct the whole set of reduced density matrices 
at weak coupling (see 24]): it turns out that the off-diagonal elements of these 
matrices decay to zero at infinity, with the same decay exponents as the non- 
interacting system; in this sense, the construction presented below rigorously 
exclude the presence of long range order in the ground state, and the absence 
of anomalous critical exponents (in other words, the interacting system is in 
the same universality class as the non- interacting one). Let me also mention 
that more sophisticated extensions of the methods exposed here also allowed 
us to: (i) give an order by order construction of the ground state correlations 



26, 



of the same model in the presence of electromagnetic interactions 
predict a possible mechanism for the spontaneous generation of the Peierls-Kekule 
instability [27|; (iii) prove the universality of the optical conductivity (28|. 
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The plan of these lectures is the following: 



• In Section |2l I introduce the model and state the main result. 



• In Section |3l I review the non-interacting case. 



• In Section HJ I describe the formal series expansion for the ground state 
energy, I explain how to conveniently re-express it in terms of Grassmann 
functional integrals and estimate by naive power-counting the generic A^-th 
order in perturbation theory, so identifying two main issues in the conver- 
gence of the series: a combinatorial problem, related to the large number 
of Feynman graphs contributing at a generic perturbative order, and a di- 
vergence problem, related to the (very mild) ultraviolet singularity of the 
propagator and to its (more serious) infrared singularity. 



• In Section |5l I describe a way to reorganize and estimate the perturbation 
theory, via the so-called determinant expansion, that allows one to prove 
convergence of the series, for any fixed choice of the ultraviolet and infrared 
cut-offs. 



In Section [6l I describe how to resum the determinant expansion in order 
to cure the mild ultraviolet divergences appearing in perturbation theory, 
via a multiscale expansion, whose result is conveniently expressed in terms 
of Gallavotti-Nicold trees 20, 21 



• In Section [TJ I describe how to resum and cure the infrared divergences and 
conclude the proof of the main result. 



• Finally, in Section |8l I draw the conclusions. A few technical aspects of the 
proof are described in the Appendixes 



The material presented in this lectures is mostly taken from 2J]. Some tech- 
nical proofs concerning the determinant and the tree expansions are taken from 
the review 23[. Other reviews of the constructive RG methods discussed here, 
which the reader may find useful to consult, are 

i3,S,|39|,|43|- 
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2. THE MODEL AND THE MAIN RESULTS 

The grandcanonical Hamiltonian of the 2D Hubbard model on the honeycomb 
lattice at half filling in second quantized form is given by: 

j=l,2,3 

where: 

1. Aa = a is a periodic triangular lattice, defined as A = B/LB, where L G 
N and B is the infinite triangular lattice with basis h = |(3, -\/3), h = 
I (3, —y/S). Ab = Aa+Si is obtained by translating A^ by a nearest neighbor 
vector 6i, i = 1,2, 3, where 

5i = (l,0), 62 = ^{-l,V3), 53 = i(-l,-V3). P2) 

The honeycomb lattice we are interested in is the union of the two triangular 
sublattices A^ and A^, see FigJT] 

2. ai^ are creation or annihilation fermionic operators with spin index a =^4- 
and site index x G A^, satisfying periodic boundary conditions in x. Sim- 
ilarly, fei^ are creation or annihilation fermionic operators with spin index 
a =ti and site index x G A^, satisfying periodic boundary conditions in x. 

3. U is the strength of the on-site density-density interaction; it can be either 
positive or negative. 

Note that the terms in the second line of Eq. (l2[ip can be rewritten as the sum of 
a truly quartic term in the creation/annihilation operators (the density- density 
interaction), plus a quadratic term (a chemical potential term, of the form —fiN, 
with /i = f//2 the chemical potential and the total particles number operator), 
plus a constant (which plays no role in the thermodynamic properties of the 
system). The Hamiltonian fl2[ll) is hole-particle symmetric, i.e., it is invariant 
under the exchange < — >ai , < — > — 6^ , . This invariance implies 
in particular that, if we define the average density of the system to be p = 
(2|A|)-^(A^). with = (at aZ +bt r ? ) the total particle number 
operator and (■)^^ = Tr{e~^^^-}/Tr{e~''^'^} the average with respect to the 
(grandcanonical) Gibbs measure at inverse temperature fi, one has p = 1, for 



FIG. 1. A portion of the honeycomb lattice A. The white and black dots correspond to 
the sites of the A^ and A^ triangular sublattices, respectively. These two sublattices 
are one the translate of the other. They are connected by nearest neighbor vectors 
^ii ^3 that, in our units, are of unit length. 

any |A| and any /3; in other words, the grand-canonical Hamiltonian Eg. (1211 p 
describes the system at half-filling, for all U,P,A. Let 

MU) = ~ \im |A|-MogTr{e-^^n- W 

p |A|-s>oo 

be the specific free energy of the system and e{U) = lim^_>.oo //3(?7) the specific 
ground state energy. We will prove the following Theorem. 

Theorem 2.1 There exists a constant Uq > such that, if \U\ < Uq, the specific 
free energy fp{U) of the 2D Hubbard model on the honeycomb lattice at half filling 
is an analytic function of U, uniformly in /3 as (3 ^ oo, and so is the specific 
ground state energy e{U). 

The proof is based on RG methods, which will be reviewed below. A straight- 
forward extension of the proof of Theorem 12.11 allows one to prove that the cor- 
relation functions (i.e., the off-diagonal elements of the reduced density matrices 
of the system) are analytic functions of U and they decay to zero at infinity with 



the same decay exponents as in the non- interacting {U = 0) case, see [24] • This 
rigorously excludes the presence of LRO in the ground state and proves that the 
interacting system is in the same universality class as the non-interacting system. 
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3. THE NON-INTERACTING SYSTEM 

Let us begin by reviewing the construction of the finite and zero temperature 
states for the non-interacting [U = 0) case. In this case the Hamiltonian of 
interest reduces to 

1=1,2,3 



with A, ai^, ? defined as in items (l)-(4) after (I2[ip . We aim at computing 
the spectrum of by diagonahzing the right hand side (r.h.s.) of f l3[ip . To this 
purpose, we pass to Fourier space. We identify the periodic triangular lattice A 
with the set of vectors of the infinite triangular lattice within a "box" of size L, 
i.e., 

A = Kfi + riafa : < rii, ^2 < L - 1} , §^2) 

with li = |(3, v^) and I2 = |(3, — v^). The reciprocal lattice A* is the set of 
vectors K such that e*^^ = 1, if x G A. A basis Gi, G2 for A* can be obtained by 
the inversion formula: 

'^ii Gi2\ ^ 2^ All 
G21 G22 J y-u I22 j 

which gives 

Gi = y(l,v^), G2 = y(l,-v^). P4) 
We call Bl the set of quasi-momenta k of the form 

^=^Gi + ^G2, mi,m2GZ, PS) 

identified modulo A*; this means that Bl can be identified with the vectors k of 
the form (121 2p and restricted to the first Brillouin zone: 

Bl = {k='^Gi + '^G2 : 0<mi,m2<L-l} . §^6) 

Given a periodic function / : A — )■ M, its Fourier transform is defined as 



which can be inverted into 

f{k) = j2e-'''fix), keBL, 
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where we used the identity 

and 6 is the periodic Kronecker delta function over A*. 



corresponding set of operators in momentum space: 



We now associate to the set of creation/annihilation operators ai^, b^_^g ^ the 



a 

|A| 



Using dUD-dlS]), we find that 



:,=E»"*<'t, ''t=E»"*^+f,. 011) 

are fermionic creation/annihilation operators, periodic over A*, satisfying 
and {a| , 6g ,} = 0- With the previous definitions, we can rewrite 

= E E (4.&,,4. + &^,4.%,J = pis) 

xeA (T=t4- 

i=l,2,3 



|A| A l^^(A^) iU 



fc,(T 



<7=n 



A:,. 



with vo = |t the unperturbed Fermi velocity (if the hopping strength t is chosen 
to be the one measured in real graphene, Vq turns out to be approximately 300 
times smaller than the speed of light) and 



n{k) = - ^ e'^^^'^'^^ = - 1 + 2e-'i'=i cos( — fcs) (12114) 



3 ^ 3 



V3, 



=1 



2 



the complex dispersion relation. By looking at the second line of Eq.( ]3[13ll . 
one realizes that the natural creation/annihilation operator is a 2D spinor of 
components a and b: 



whose real space counterparts read 



*t = («t.ft«,j. n,= [i^': I ■ 016) 
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In order to fully diagonalize the theory, one needs to perform the diagonalization 
of the 2x2 quadratic form in the second line of Eq. (131 13p , which can be realized 
by the fc-dependent unitary transformation 

in terms of which, defining 



V2|!^(A:)| V2 



^^k,a/ y \n{k)\ k,a^ k,. 



we can rewrite 



A A ^ k,a k k I Q } k k k,a 

k&BL ^ 

<T=n 



|A| '^.^ \ -\Vl{k)\ 

keBL 



keBL 

<T=n 

The two energy bands ±fo|r2(/c)| are plotted in FigJJl They cross the Fermi 
energy Ep = at the Fermi points k = pp, u = ±, with 

Pf = (y,^^), P20) 

close to which the complex dispersion relation vanishes linearly: 

Q{p^ + k') = ik[ + uk'^ + 0{\k'\^) , ([3121) 

resembling in this sense the relativistic dispersion relation of (2 + l)-dimensional 
Dirac fermions. From Eq. fl3[19p it is apparent that the ground state of the system 
consists of a Fermi sea such that all the negative energy states (the "a-states" ) are 
filled and all the positive energy states (the "/3-states") are empty. The specific 
ground state energy eo,A is 



keBL 



from which we find that the specific ground state energy in the thermodynamic 
limit is 

r flk 

e{0) = -2vo / jj^Mk)\, (1123) 
Jb \^\ 





FIG. 2. A sketch of the energy bands of the free electron gas with nearest neighbor 
hopping on the honeycomb lattice. The red plane corresponds to the Fermi energy at 
half-filling. It cuts the bands at a discrete set of points, known as the Fermi points 
or Dirac points. From the picture, it seems that there are six distinct Fermi points. 
However, after identification of the points modulo vectors of the reciprocal lattice, it 
turns out that only two of them are independent. 

where B := {k = hGi + ^2^2 : t, e [0, 1)} and \B\ = 8717(3^3). Similarly, the 
finite volume specific free energy /q^ := — (/3|A|)~^ logTrje"^^*} is 



0,A 



P\A\ 



^log[(l + e^''ol^(^^)l)(l + 



-/3t,oif^(fc)n 



(E124) 



from which 



MO) 



2 r dk ^ 



(1+ 



,l3vo\n{k)\ 



)(i + 



-l3vo\n{k)\ 



Besides these thermodynamic functions, it is also useful to compute the Schwinger 
functions of the free gas, in terms of which, in the next section, we will write down 
the perturbation theory for the interacting system. These are defined as follows. 
Let xo G [0,/3) be an imaginary time, let x := (xo,x) G [0,(3) x A and let 
us consider the time-evolved operator \l/^o- = e^^°\l/i^e~^^°, where is the 
two-components spinor defined in Eq. (l3[16p : in the following we shall denote its 
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components by ^^J^^, with p e {1,2}, -^^^ ^ = and ^^^^,2 = K+5^,a (here 
^1 = (0,(5i)). We define the n-points Schwinger functions at finite volme and 
finite temperature as: 

where: Xj G [0, f3] x A, cr^ =ti, = ±, Pi = 1, 2 and T is the operator of fermionic 
time ordering, acting on a product of fermionic fields as: 

where vr is a permutation of {l,...,n}, chosen in such a way that X7r(i)o > 
■ ■ ■ > a;^(„)o, and (—1)'^ is its sign. [If some of the time coordinates are equal each 
other, the arbitrariness of the definition is solved by ordering each set of operators 
with the same time coordinate so that creation operators precede the annihila- 
tion operators.] Taking the limit A — )■ oo in (131 26p we get the finite tempera- 
ture n-point Schwinger functions, denoted by 5'^(xi, ei, cxi, pi; . . . ; x„. En, Cn, Pn), 
which describe the properties of the infinite volume system at finite temper- 
ature. Taking the /3 — )■ cxd limit of the finite temperature Schwinger func- 
tions, we get the zero temperature Schwinger functions, simply denoted by 
5'„(xi, £1, (Ti, pi; . . . ; x„, £„, cr„, p„), which by definition characterize the proper- 
ties of the thermal ground state of fl2[ll) in the thermodynamic limit. 

In the non-interacting case, i.e., if H/^ = H^, the Hamiltonian is quadratic in 
the creation/annihilation operators. Therefore, the 2n-point Schwinger functions 
satisfy the Wick rule, i.e., 

G^.=^^.'^^(T{^x-,....*,^,,.^,J)^^^. il28) 

Moreover, every n-point Schwinger function S'^'^(xi, ei, ai, pi; . . . ; x„, £„, cr„, p„) 
with ^27=1 7^ is identically zero. Therefore, in order to construct the whole 
set of Schwinger functions of H^, it is enough to compute the 2-point function 
S'q'^(x — y) = {T^{'^^^a,p^y a p'}) j3 \- '^^^^ '^^^ casily reconstructed from the 
2-point function of the a-fields and /3- fields, see Eq. (131 181) . 

Let X e A, = lAI"^ Efes, and = |A|-i Efce. if 

X = (xo,x), let a^^^ = e^A^'oaJ^e--^A^o and l3^^^ = e^^^o p±^f.-K^o ^ \ straight- 
forward computation shows that, if —f3 < Xq — yo < (3, 

kdBL 

P ^VQ{xQ-yo)\Q.{k)\ ^vo(xQ-yo+l3)\Q.(k)\^ 

• l(a;o -z/0 > 0)- — — - U(xo -1/0 < 0)— ^ , 
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1 (^0 - 1/0 > 0) ^ , - a (^0 - 1/0 < 0) 



and (T{a,,,/3+^,})^_^ = (T{/3,-,«+^,})^_^ = 0. A priori, Eqs.dSISni) and (EE 
are defined only for —(3 < Xq — yo < (3, but we can extend them periodically 
over the whole real axis; the periodic extension of the propagator is continuous in 
the time variable for xq — yo ^ and it has jump discontinuities at the points 
^0 " Vo ^ Note that at Xq — y^ = f3n, the difference between the right and 
left limits is equal to (— l)"fe,y-, so that the propagator is discontinuous only at 
X - y = /3Z X 0. If we define B^^l ■= Bp x Bl, with Bp = {ko = ^{no + ^ : 
no G Z}, then for x — y ^ /3Z x we can write 

(T{«-x.'}),A = m\ 5Z ^-^"^^"^^ , ^ ' (1331) 

^ ^l^^'k^.i -«A;o + t;o|^^(A^)| 

If we now re-express a^^^ and /3x,o- in terms of a^^^ and &x+5i o-' using Eq.f l3[18p . 
we find that, for x — y ^ /3Z x 0: 



So^i^ - y)p,p' ■■= S!^'^{x, a, -, p; y, a, +, p') 



(7=0 



Finally, if x - y = (0^, 0): 



4. PERTURBATION THEORY AND GRASSMANN INTEGRATION 



Let us now turn to the interacting case. The first step is to derive a formal 
perturbation theory for the specific free energy and ground state energy. In other 
words, we want to find rules to compute the generic perturbative order in U of 
fp^^ ■= _(/3|A|)-i logTr{e-^-^A}_ We write i^A = + '^a, with Vk the operator 
in the second line of Eq. (I2[ip and we use Trotter's product formula 

(HI) 



-'^^^ = hm e-^^A/"(l - ^Vj 

n— >oo . n, 
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FIG. 3. The four- legged graph element (left); an example of a Feynman diagram of 
order 3 (right). 



so that, defining V\{t) := e aV^c 



A 



Tr{e-^^A} _ 
Tr{e-^^S} ~ 

/•*! /-^iv-i Trlp-I^K 



(!12) 



iV>l 



Using the fermionic time-ordering operator defined in Eg. (131271) . we can rewrite 
Eq.dlE]) as 

where (■)° = Tr{e-'3^^X.}/Tr{e-'^^X}, 

and Jj.^ (ix must be interpreted as Jj.^ (ix = /_^/2 '^^o S^eA- Note that the 
A^-th term in the sum in the r.h.s. of Eq. (l4[3p can be computed by using the 
Wick rule (I3[28p and the explicit expression for the 2-point function Eqs. (l3[33|) - 
(131 34p . It is straightforward to check that the "Feynman rules" needed to compute 
^ are the following: (i) draw N graph elements consisting of 
4- legged vertices, with the vertex associated to two labels Xj and pi, i = 1, . . . , N, 
and the four legs associated to two exiting fields (with labels (xj, f, pi) and (x,, J, 
, Pi)) and two entering fields (with labels (xj,'|', p^) and (xj,J,, p^)), respectively; 
(ii) pair the fields in all possible ways, in such a way that every pair consists of one 
entering and one exiting field, with the same spin index, see FigJHl (iii) associate 
to every pairing a sign, corresponding to the sign of the permutation needed to 
bring every pair of contracted fields next to each other; (iv) associate to every 
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paired pair of fields [^x,^(Ti,pi; ^x,,crj ,pj oriented line connecting the z-th with 
the j-th vertex, with orientation from j to i; (v) associate to every oriented line 
[j — !■ i] a value equal to 

(!15) 

where = ^^'^ x Bl, = B^ n {ko : Xo(2-*'|A;o|) > 0} and Xo{t) is 

a smooth compact support function that is equal to 1 for (say) \t\ < 1/3 and 
equal to for \t\ > 2/3; (vi) associate to every pairing (i.e., to every Feyn- 
man graph) a value, equal to the product of the sign of the pairing times f/^ 
times the product of the values of all the oriented lines; (vii) integrate over Xj 
and sum over pi the value of each pairing, then sum over all pairings; (viii) 
finally, take the M — )• oo limit: the result is equal to (T{(V}3^a(\E'))^})° ^. Note 
that the M — i- oo limit of the propagator g{x) is equal to S'o'^(x) if x 7^ 0, 
while limM^oo5'(0) = S'o'^(0~,0) + ^, see Eq.f l3l34p : the difference between 
limM->oo fi'(x) and Sq'^{'x) takes into account the — | terms in the definition of 

An algebraically convenient way to re-express Eq.f l4[3p is in terms of Grass- 
mann integrals. Consider the set Am,i3,l = {'^kcrp}'^"^!^^"^'^' where the Grass- 

mann variables ip^^p satisfy by the definition the anticommutation rules {"^Ako-p' 
'^k',a',p'} — 0- particular, the square of a Grassmann variable is zero and the 
only non-trivial Grassmann monomials are at most linear in each variable. Let 
the Grassmann algebra generated by Am,i3,l be the set of all polynomials ob- 
tained by linear combinations of such non-trivial monomials. Let us also define 
the Grassmann integration J [HkeB'"' 11^=^ '^^ko-p'^^ko-p] linear oper- 

ator on the Grassmann algebra such that, given a monomial Q{'ip~ ^-ij)^) in the 
variables ^k,(T,p5 action on Q{'ip~ ^-ip^) is except in the case Q{^~,-4'~^) = 
HkeB^") nCT=t4 '^kcrp'V^ko-p; ^ permutation of the variables. In this case the 

value of the integral is determined, by using the anticommutation properties of 
the variables, by the condition 

p=l,2 p=l,2 

/ [ n n di^t^AK<r,p\ n n ^k^^jlp = ^ 

Defining the free propagator matrix ^k as 



—Vo\l{k) —tkQ 
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and the "Gaussian integration" PMidip) as 



/3,-L 



exp { - (/3|Ari ^t^9^'^K^} ' 



(A/) 
'/3,-L 



it turns out that 



/ Pidi^)^^,,^M,^. = /3|A|<5.„..4„k.^k, , (119) 

while the average of an arbitrary monomial in the Grassmann variables with 
respect to PM{dip) is given by the fermionic Wick rule with propagator equal to 
the r.h.s. of Eq. fl4[9p . Using these definitions and the Feynman rules described 
above, we can rewrite Eq. fl4[3p as 



J I = lim /pM(#)e-^('^) , illO) 
Tr|e~^ a| m-s>oo J 

where 



p=l,2 • 

^t,P = ^ E ^^''^^t,P' xG(-/3/2,/3/2]xA (1112) 



and the exponential e"^*-'''^ in the r.h.s. of Eq. fl4[10p must be identified with its 
Taylor series in U (which is finite for every finite M, due to the anticommutation 
rules of the Grassmann variables and the fact that the Grassmann algebra is finite 
for every finite M). Apriori, Eq. fl4[10p must be understood as an equality between 
formal power series in U. However, it can be given a non-perturbative meaning, 
provided that we can prove the convergence of the Grassmann functional integral 
in the r.h.s., as shown by the following Proposition. 

Proposition 4.1 Let 

■■= log / PMidi^) (e-^W) (HIS) 

and let (5 and |A| he sufficiently large. Assume that there exists f/o > such that 
Pjl^f^ is analytic in the complex domain \U\ < Uq and is uniformly convergent as 
M -> oo. Then, z/|f/| < Uq, 

fp,A = E (2 + 2cosh{/3vom)\)) + Jim^F^^y • 0114) 
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Proof. We need to prove that 



rr^L^nJ^ = exp { - /3|A| Jim F^^ } (ill5) 



under the given analyticity assumptions on -Fg^ • The first key remark is that, 
if /3,A are finite, the left hand side of Eq.f l4[15|) is a priori well defined and 
analytic on the whole complex plane. In fact, by the Pauli principle, the Fock 
space generated by the fermion operators a^^, 6^^- ^, with a; G A, cr =ti, is finite 
dimensional. Therefore, writing H\ = + Va, with and Va two bounded 
operators, we see that Tr{e~^^'^} is an entire function of U, simply because e~^^^ 
converges in norm over the whole complex plane: 



I \n-k 



°" on °° /^fc I It/ I |fc flri—k\ I TT 

n=0 k=0 n>k ^ ' 

= e^ll^All+/^II^A|| ^ (15116) 

where the norm || ■ || is, e.g., the Hilbert-Schmidt norm ||y4|| = ^J'Yxi^A^ A). 

On the other hand, by assumption, F^^^ is analytic in |f/| < f/g, with Uq inde- 
pendent of /3, A, M, and uniformly convergent as M — )■ oo. Hence, by Weierstrass 
theorem, the limit F/j^a = liniAf^oo -F^^'' is analytic in |f/| < Uq and its Taylor 
coefficients coincide with the limits as M — )■ oo of the Taylor coefficients of Fj^\ 

Moreover, limAf-s>oo e"^'^'^"-'^ = e"''''^'"^''-^, again by Weierstrass theorem. 

As discussed above, the Taylor coefficients of e~^\^\^PA coincide with the Tay- 
lor coefficients of Tr{e-^^'^}/Tr{e-^-^o.'^}: therefore, Tr{e-^-^A}/Tr{e-^-^0'A} = 
g-/3|A|F^,A the complex region |f/| < f/g, simply because the l.h.s. is entire in 
f/, the r.h.s. is analytic in |f/| < [/q and the Taylor coefficients at the origin of 
the two sides are the same. Taking logarithms at both sides proves Eg. ( ]4[14l) . ■ 

By Proposition 14.11 the Grassmann integral Eg. (141 131) can be used to com- 
pute the free energy of the original Hubbard model, provided that the r.h.s. of 
Eg.f l4[13|) is analytic in a domain that is uniform in M, /3, A and that it converges 
to a well defined analytic function uniformly as M — i- oo. The rest of these notes 
are devoted to the proof of this fact. We start from Eg.f l4[13l) . which can be 
rewritten as 

where the truncated expectation S'^ is defined as 

niV f 
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More in general, 



Q7V r 

and £^"^(Vi, . . . , Vat) L „ = £^{V;N). It can be checked by induction that 



iv,=v 

the truncated expectation is related to the simple expectation £{X{ip)) = 
J PM{d^)X{ij) by 

N 

m=i (yi,...,y'") 

where the second sum in the r.h.s. runs over partitions of {1, . . . , A^} of multi- 
plicity m, i.e., over m-ples of disjoint sets such that U™]^F* = {1, . . . ,N}, with 
= {jl, . . . ,j\Yi\}- Note that £{V'^) = SiVi, . . . , Vat)]^^^ can be computed as 
a sum of Feynman diagrams whose values are determined by the same Feynman 
rules described after Eq. fl4[4l) (with the exception of rule (viii): of course, since 
S{X) = J PM{dip)X, M should be temporarily kept fixed in the computation); 
we shall write 

where Fat is the set of all Feynman diagrams with vertices, constructed with 
the rules described above; Val(^) includes the integration over the space-time 
labels Xj and the sum over the component labels pi: if ^ G F^, we shall write 
symbolically 

Val(^) = agf/^ J2 [ d^i---d^^l[5^^e),a'ii)9p(e),p'(eM()-^'ii)) , (!122) 

P1,...,PN ^65 

where ag is the sign of the permutation associated to the graph Q and we denoted 
by (x(£), cr(£), p(£)) and (x'(£), cr'(£), p'(£)) the labels of the two vertices, which 
the line £ exits from and enters in, respectively. Using Eqs. fl4[20p - fl4[2ip . it can 
be proved by induction that 

S^{V;N)= 5^ Val(e;), (1123) 

where F^ C F^r is the set of connected Feynman diagrams with N vertices. Com- 
bining Eg. (141 171) with Eq. (l4[23l) we finally have a formal power series expansion 
for the specific free energy of our model (more precisely, of its ultraviolet regu- 
larization associated to the imaginary-time ultraviolet cutoff xo{'2~^^\ko\))- The 
Feynman rules for computing Val(^) allow us to derive a first very naive upper 
bound on the A^-th order contribution to F^^[^ , that is to 

C^"^--^n7r^"(V;A^)- il24) 
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We have: 

I3\A.\N\ 



<^2^|f;|A.||^||^+i||^||A.-i^ (1225) 

where |r^| is the number of connected Feynman diagrams of order and 
/3|A|(2|[/|)^||5f||^+^ ll^fll^""'^ is a uniform bound on the value of a generic con- 
nected Feynman diagram of order A^. The bound is obtained as follows: given 
Q G r^, select an arbitrary "spanning tree" in i.e. a loopless subset of Q 
that connects all the N vertices; now: the integrals over the space-time coordi- 
nates of the product of the propagators on the spanning tree can be bounded 
by /3|A| llsfllf^"^; the product of the remaining propagators can be bounded by 
||5f||^+^; finally, the sum over the pi labels is bounded by 2^. Using Eq. fl4[25p 
and the facts that, for a suitable constant C > 0: (i) |r^| < C^(A^!)^ (see, e.g.. 
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Appendix A. 1.3] for a proof of this fact), (ii) \\g\\oo < CM, (iii) \\g\\i < CP, 
we find: 

IF^^/'^^^I < {2C^\U\fN\ . (1226) 

Remark. While the bound \\g\\i < C(3 (see Appendix lAl for a proof) is 
dimensionally optimal, the estimate \\g\\oo < CM could be improved to Hs^Hoo < 
(const.), at the cost of a more detailed analysis of the definition of 5'(x), which 
shows that the apparent ultraviolet logarithmic divergence associated to the sum 
over ko in Eg. (141 5 p is in fact related to a jump singularity of g{xo, 0) at xq = /3Z. 
This can be proved along the lines of j^. Section 2]. However, to the purpose of 
the present discussion, the rough (and easier) bound ||5'||oo < CM is enough; see 
Appendix El for a proof. 

The pessimistic bound Eg. (14126]) has two main problems: (i) a combinatorial 
problem, associated to the A^!, which makes the r.h.s. of Eg. ( 14126]) not summable 
over A^, not even for finite M and (3; (ii) a divergence problem, associated to 
the factor M^^^/?^"-*^, which diverges exponentially as M — )■ cxd (i.e., as the 
ultraviolet regularization is removed) and as /3 — cxd (i.e., as the temperature 
is sent to 0). The combinatorial problem is solved by a smart reorganization 
of the perturbation theory, in the form of a determinant expansion, together 
with a systematic use of the Gram-Hadamard bound. The divergence problem is 
solved by systematic resummations of the series: we will first identify the class of 
contributions that produce ultraviolet or infrared divergences and then we show 
how to inductively resum them into a redefinition of the coupling constants of the 
theory; the inductive resummations are based on a multiscale integration of the 
theory: at the end of the construction, they will allow us to express the specific 
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free energy in terms of modified Feynman diagrams, whose values are not affected 
anymore by ultraviolet or infrared divergences. 

5. THE DETERMINANT EXPANSION 

Let us now show how to attack the first of two problems that arose at the 
end of previous section. In other words, let us show how to solve the combina- 
torial problem by reorganizing the perturbative expansion discussed above into 
a more compact and more convenient form. In the previous section we discussed 
a Feynman diagram representation of the truncated expectation, see Eq. fl4[23l) . 
A slightly more general version of Eq.f l4[23|) is the following. For a given set of 
indices P = {/i, . . . , /|p|}, with fi = {xi,ai, pi,ei), G {+,-}, let 

Each field a{f) p{f) ^^"^ represented as an oriented half-line, emerging from 
the point x(/) and carrying an arrow, pointing in the direction entering or exiting 
the point, depending on whether e{f) is equal to — or +, respectively; moreover, 
the half-line carries two labels, cr(/) G {tii} and p(/) G {1,2}. Now, given s 
set of indices Pi, . . . ,Ps, we can enclose the points x(/) belonging to the set Pj, 
for some j = 1, . . . , s, in a box: in this way, assuming that all the points x(/), 
/ G UjPj, are distinct, we obtain s disjoint boxes. Given V := {Pi, . . . , Ps}, we 
can associate to it the set r^('P) of connected Feynman diagrams, obtained by 
pairing the half-lines with consistent orientations, in such a way that the two 
half-lines of any connected pairs carry the same spin index, and in such a way 
that all the boxes are connected. Using a notation similar to Eq. (141 221) . we have: 

Val(6;) = ^^11 ^<r(i)yie)9p(e),p'ii)W) - x'(£)) , (02) 

A different a more compact representation for the truncated expectation, alter- 
native to Eq. fl5[2p . is the following: 

£^(V'P,,...,V^pJ= Yl (^rllge f dPT{t)detG^{t) , (ElS) 

where: 

• any element T of the set T{V) is a set of lines forming an anchored tree 
between the boxes Pi, . . . , Ps, i.e., T is a set of lines that becomes a tree if 
one identifies all the points in the same clusters; 
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• a-r is a sign (irrelevant for the subsequent bounds); 

• gi is a shorthand for 5a{e),cT'ie) gp(e),p'{e)W) - 

• if t = {ti^ii G [0, 1], 1 < i, i' < s}, then dPT^t) is a probabihty measure with 
support on a set of t such that tj^j/ = Uj ■ Uj/ for some family of vectors 
Uj G M" of unit norm; 

• if 2n = Yli=i then G"^(t) is a (n — s + 1) x (n — s + 1) matrix, whose el- 
ements are given by G^^, = ti(f),i{f')ge{fj'), where: /, /' ^ Ui^xifl, fj] and 

fj are the two field labels associated to the two (entering and exiting) 
half-lines contracted into £; i{f) G {1, . . . , s} is s.t. / G Pi{f)', 9iifj') is the 
propagator associated to the line obtained by contracting the two half-lines 
with indices / and /'. 

If s = 1 the sum over T is empty, but we can still use the Eq. fl5[3p by in- 
terpreting the r.h.s. as equal to 1 if Pi is empty and equal to det G^{1) otherwise. 



The proof of the determinant representation is described in Appendix [Bl this 
representation is due to a fermionic reinterpretation of the interpolation formu- 
las by Battle, Brydges and Federbush jl, 10, [ll|, used originally by Gawedski- 
Kupianen 22| and by Lesniewski js^, among others, to study certain (1 + 1)- 
dimensional fermionic Quantum Field Theories. Using Eq.f l5[3p we get an alter- 
native representation for the A^-th order contribution to the specific free energy: 



,(M;Af) 



-1) 



N 



(V; N) 



■1) 



N 



Pi,...,PN TeTff 



dxi • • • dxjv 



/3\A\ m ' ' ' /3\A\ m 

■n^-W--'W^pW.p'w(^(^) -^'W) / ^^r(t)detG^(t) . ^4) 



Using the fact that the number of anchored trees in T^r is bounded by C^Nl 
for a suitable constant C (see, e.g., (23|, Appendix A. 3. 3] for a proof of this fact), 
from Eq. fl5[4p we get: 



r^{M;N) I 
"l3,A I 



< {const. y \U\'' ||(7||f-^||detG 



(05) 



In order to bound detG-^, we use the Gram-Hadamard inequality, stating that, if 
M is a square matrix with elements Mij of the form Mij = {Ai,Bj), where Ai, 
Bj are vectors in a Hilbert space with scalar product (■, ■), then 



detM| < Yl\\Ai\ \ ■ \\Bi 



(06) 



where || ■ || is the norm induced by the scalar product. See 23|, Theorem A.l] for 
a proof of Eq. (l5[6p . 
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Let = ® T-Lq, where "Hq is the Hilbert space of the functions F : 
[-(3/2,(3/2) X A ^ C^ with scalar product (F, G) = Ep=i,2 / dz F*{z)Gp{z), 
where Fp = [F]p, Gp = [G]p, p = 1,2, are the components of the vectors F and 
G. It is easy to verify that 

^/j' = Mf),i{f') ^TW.a'w gp[f),p(f){y^U) - = 

where: Uj G M", i = l,...,n, are vectors such that tj^j/ = Uj ■ Uj/; e-f- = (1,0), 
Bi = (0, 1); Ax,p and Bx,p have components: 



1 ^ v/xo(2-^|A;o|)e-'''^(^-'') f tk^ -voQ*{k) 



so that 



keB}j^ L-^u ' -ui--\-/i J \ / i,p 



^' ' kPB(*^' ^0 + ^ol^(^)l 

for a suitable constant C. Using the Gram-Hadamard inequality, we find 
1 1 det G"^! loo < (const. )^M^+^; substituting this result into Eq. fl5[5l) . we finally 
get: 

li^jX^^^I < (constO^I^/rM^+V"^-' , (010) 

which is similar to Eq. (141261) . but for the fact that there is no A^! in the r.h.s.! 
In other words, using the determinant expansion, we recovered the same dimen- 
sional estimate as the one obtained by the Feynman diagram expansion and we 
combinatorially gained a 1/A^!. The r.h.s. of Eq. flSllOp is now summable over N 
for \U\ sufficiently small, even though non uniformly in M and /3. In the next 
section we will discuss how to systematically improve the dimensional bound by 
an iterative resummation method. 



6. THE MULTISCALE INTEGRATION: THE ULTRAVIOLET 

REGIME 

In this section we begin to illustrate the multiscale integration of the fermionic 
functional integral of interest. This method will later allow us to perform iterative 
resummations and to re-express the specific free energy in terms of a modified 
expansion, whose A^-th order term is summable in and uniformly convergent 
as M — )■ oo and (3 — oo, as desired. 
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The first step in the computation of the partition function 



and of its logarithm is the integration of the ultraviolet degrees of freedom corre- 
sponding to the large values of k^. We proceed in the following way. We decom- 
pose the free propagator into a sum of two propagators supported in the regions 
of ko "large" and "small" , respectively. The regions of ko large and small are de- 
fined in terms of the smooth support function Xo(^) introduced after Eq. fl4[5l) : 

note that, by the very definition of Xo; the supports of Xo(y\J k^ + \k — j and 

XoiysJ kl + \k — j)]p\'^ j are disjoint (here | ■ | is the euchdean norm over M^/A*). 
We define 



A.,.(k) = 1 -Xo(V^o' + \k-i^F? )-X^{\Jkl + {k-p-p? ) (02) 
and /j..r.(k) = 1 — fu.v.i}^)-, so that we can rewrite as: 

g^. = /n...(k)^k + A..(k)^k =^^("-''-)(k) + g^'-^-\k) . ^3) 



We now introduce two independent set of Grassmann fields {ipi^^p"^} and 
{'^ka^ft}' with k G B^^2 > ^ =ti) P = 1;2, and the Gaussian integrations 
P(d^'("-''-)) and Pidip^'-''-^) defined by 



I p(rf^('-'--))C;!"v^£:^ = ws..,.,s^.MJ-\k,) . (114) 

Similarly to PMidijj), the Gaussian integrations P{dip^'^''"'^), P{d-ip^^''^'^) also admit 
an explicit representation analogous to (I4[8p . with ^k replaced by ^{«.^.)(k) or 
^(*-'--)(^k) and the sum over k restricted to the values in the support of fu.v.(k) 
/j.r.(k), respectively. The definition of Grassmann integration implies the follow- 
ing identity ("addition principle"): 

P(#)e-^(^) = [ P(#(*-'--)) [ p(rf^(--))e-v(^(-'-)+^(--)) 



so that we can rewrite the partition function as 

-/3|A|f/*^) 



■ ,9, A 



f P(#(--))exp {J2^_^u.A-V{^^'-'-^ + -y,n)} 

n>l 



p(rf^/'(""^)e-^«^'^" , die) 
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where the truncated expectation £j„ is defined, given any polynomial Vl{'^p^'^^'"^^) 
with coefficients depending on ip^'^-^-\ as 



A=0 



and Vo is fixed by the condition Vo(0) = 0. We will prove below that Vq can be 
written as 

oo n 

v„w = E E E [IKi.^-.CU; ■ 

n=l o-i,...,o-„=ti Pl,.-->P2n=l,2 kl,...,k2n i=l 

n 

■WM,2nAkl, k2n-l) ^(^^(ks,-! " ka,)) , ^8) 

where p = (pi, . . . , p2n) and we used the notation 
6{k) = 6{k)6ik,) , 6{k) = \A\ kn.G^+n.G,^ S{ko) = P6,,,o , (09) 

ni,n2ez 

with Gi,G2 a basis of A*. The possibility of representing Vm in the form (I6[8I) . 
with the kernels WM.2n,p independent of the spin indices (Tj, follows from a number 
of remarkable symmetries, discussed in Appendix [Cl see in particular symmetries 
(l)-(3) in Lemma [C.ll The regularity properties of the kernels are summarized 
in the following Lemma, which will be proved below. 

Lemma Ol The constant Fq^m in Ii6[6\) and the kernels WM,2n,p in Ii6\8\) are 
given by power series in U, convergent in the complex disc \U\ < Uq, for Uq 
small enough and independent of M, /3, A; after Fourier transform, the ^-space 
counterparts of the kernels WM,2n,p satisfy the following bounds: 

j dXi---dX2„|l^M,2n,p(Xl,...,X20| </3|A|C'"|t/r'^^^l'"-l> (EllO) 

for some constant C > 0. Moreover, the limits Fq = limjv/-5>oo -^o,m o,nd 
W^2n,p(xi, . . . , X2n) = limAf^oo WAf,2n,p(xi, . . . , X2„) cxist and arc reached uni- 
formly in M, so that, in particular, the limiting functions are analytic in the 
same domain \U\ < Uq and so are their /3, |A| — )■ oo limits (that, with some abuse 
of notation, we shall denote by the same symbols). 

Remark. Once that the ultraviolet degrees of freedom have been integrated 
out, the remaining infrared problem (i.e., the computation of the Grassmann 
integral in the second line of Eq.f l6l6p ) is essentially independent of M, given the 
fact that the limit W2n,p of the kernles WM,2n,p is reached uniformly and that the 
limiting kernels are analytic and satisfy the same bounds as Eq. (l6[10p . For this 
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reason, in the infrared integration described in the next two sections, M will not 
play any essential role and, whenever possible, we will simplify the notation by 
dropping the label M. 

Before we present the proof of Lemma Hi, let us note that the kernels WM,2n,p 
satisfy a number of non-trivial invariance properties. We will be particularly 
interested in the invariance properties of the quadratic part WA,f,2,(pi,p2)(^^i; '^s), 
which will be used below to show that the structure of the quadratic part of the 
new effective interaction has the same symmetries as the free integration. The 
crucial properties that we will need are summarized in the following Lemma, 
which is proved in Appendix [Ul 

Lemma [6]. 2. The kernel 14^2 (k) := W0w,2(k, k), thought as a2 x2 matrix with 
components PV'M,2,(i,j)(k, k), satisfies the following symmetry properties: 

= (T,W2{{ko,-k,,k2))a^ = W^{{ko,-k)) = -a,W2{{-ko,k))a3 , (011) 

where ai,a2,cr3 are the standard Pauli matrices. In particular, if W^(k) = 
lim^jA|-5-oo W^2(k), in the vicinity of the Fermi point p"^ = (0,p^), with w = ±, 

iV(k) = - {if, I, + 0(|k - pi;n , 6112) 

for some real constants Zq,6o. 

Note that Eq. fl6[12p can be read by saying that, in the zero temperature 
and thermodynamic limits, the two-legged kernel has the same structure as the 
inverse of the free covariance, 5(7^ (k), modulo higher order terms in k — p"^. 
This fact will be used in the next section to define a dressed infrared propagator 
[S'Q"^(k) -|- W^(k)]~^, with the same infrared singularity structure as the free one. 
We will come back to this point in more detail. For the moment, let us turn to 
the proof of Lemma [6l 1 , which illustrates the main RG strategy that will be also 
used below, in the more difficult infrared integration. 

Proof of Lemma [6ll. Let us rewrite the Fourier transform of ^'^''■^■^(k) as 

M 

<7("-)(x) = 5^(7W(x), (1113) 

h=l 

where 
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with H^{ko) = Xo{'2-'\ko\) and, if h > 2, H^^k^) = Xoi'^-^'M) - X^i'^-^^^M) ■ 
Note that [g'^^\Q)]pp = 0, p = 1,2, and, for any integer -ft' > 0, g''^\^) satisfies 
the bound 



1 + (2'>|xo|/3 + |x|a)^ 



\9^'\^)\\ < . , .oM^^i'' , 1^1 ^ (015) 



where |xo|/3 = min„gz |xo + is the distance over the one-dimensional torus of 
length (3 and |x|a = min^^jg is the distance over the periodic lattice A (here 

B is the triangular lattice defined after Eq.f l2[ip ): see Appendix lA] for a proof. 
Moreover, g^^\x.) admits a Gram representation that, in notation analogous to 
Eq.f l5[7p . reads 

^7S(^-y) = (A5C')' dais) 

with 

[Aifl(z)], = V JfuvO<-)Hh(ko) 6pi, mi7) 

/3|A| ^^^, + vi\Q{kW]V^ [-voQik) , 

and 

m^f-^,,,[ki+vim)?Y"~ 

for a suitable constant C. Our goal is to compute 



= exp I log J P(t/V;[i'*^l)e-^('^"-^-^+'^'""')} , (019) 

where P{dijj^^'^'^^) is the fermionic "Gaussian integration" associated with the 
propagator J2hLi d^^K^) (i-e., it is the same as P{dtp^^-'"-^)); moreover, we want to 
prove that Fo,m and Vo{ip^^''^^) are uniformly convergent as M — > oo. We perform 
the integration of (1611 9p in an iterative fashion: in fact, we will inductively prove 
that for 1 < /i < M, 



where P{dip^^'^^) is the fermionic "Gaussian integration" associated with the prop- 
agator ELi^^'^Hk) and V(*^) = V; for 1 < /i< M, 

V^'^n^) = (021) 

oo „ n 

n=l p,CT j=l 
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In order to inductively prove fl6[20p - fl6[2ip we use the addition principle to rewrite 

(16122) 

where P{dip^'^^) is the fermionic Gaussian integration with propag ator ^(''^(k). 
After the integration of ip^'^^ we define 

^_V('^-i)(V>("-')+^[i'''-il)-/3|A|e;, ^ y p^^^(fe)^|^_V('^){V>("-•)+^[l.''-ll+,/.W) ^ P23) 

which proves (161201) with 

M 



Fh= 5^ . (024) 

fc=/i+i 

Let £^/f be the truncated expectation associated to P(d'ip^'^^): then we have 

Ch + V^'^-i) (^) = 5^ ^(-l)"+i^^(V('^) + ^('^)) ; n) . (025) 



n>l 



Eq. (l6[25l) can be graphically represented as in FigJH The tree in the l.h.s., con- 

h+1 




FIG. 4. The graphical representation of V*-'^ 

sisting of a single horizontal branch, connecting the left node (called the root and 
associated to the scale label h — 1) with a big black dot on scale h, represents 
V^^~^\ip)- III the r.h.s., the term with n final points represents the corresponding 
term in the r.h.s. of Eg. (16125]) : a scale label h — 1 is attached to the leftmost node 
(the root); a scale label h is attached to the central node (corresponding to the 
action of £^); a scale label h + 1 is attached to the n rightmost nodes with the 
big black dots (representing V'-''^). Iterating the graphical equation in Fig H] up 
to scale M, and representing the endpoints on scale M + 1 as simple dots (rather 
than big black dots), we end up with a graphical representation of V'-''^ in terms 
of Gallavotti-Nicold trees, see Figj5l defined in terms of the following features. 



1. Let us consider the family of all trees which can be constructed by joining 
a point r, the root, with an ordered set of n > 1 points, the endpoints of 
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FIG. 5. A tree r G TM;h,n with n = 9: the root is on scale h and the endpoints are on 
scale M + 1. 

the unlabeled tree, so that r is not a branching point, will be called the 
order of the unlabeled tree and the branching points will be called the non 
trivial vertices. The unlabeled trees are partially ordered from the root to 
the endpoints in the natural way; we shall use the symbol < to denote the 
partial order. Two unlabeled trees are identified if they can be superposed 
by a suitable continuous deformation, so that the endpoints with the same 
index coincide. It is then easy to see that the number of unlabeled trees 



with n end-points is bounded by 4"^ (see, e.g., [23!, Appendix A. 1.2] for a 
proof of this fact). We shall also consider the labelled trees (to be called 
simply trees in the following); they are defined by associating some labels 
with the unlabelled trees, as explained in the following items. 

2. We associate a label < /i < M — 1 with the root and we denote by 
TM;h,n the Corresponding set of labeled trees with n endpoints. Moreover, 
we introduce a family of vertical lines, labeled by an integer taking values in 
[h, M + 1], and we represent any tree r G TM;h,n so that, if v is an endpoint, 
it is contained in the vertical line with index = M + 1, while if it is a 
non trivial vertex, it is contained in a vertical line with index h < < M, 
to be called the scale of v; the root r is on the line with index h. In general, 
the tree will intersect the vertical lines in set of points different from the 
root, the endpoints and the branching points; these points will be called 
trivial vertices. The set of the vertices will be the union of the endpoints, 
of the trivial vertices and of the non trivial vertices; note that the root is 
not a vertex. Every vertex w of a tree will be associated to its scale label 
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hy, defined, as above, as the label of the vertical line whom v belongs to. 
Note that, if vi and V2 are two vertices and vi < V2, then h^^ < h^^. 

3. There is only one vertex immediately following the root, called vq and with 
scale label equal to h + 1. 

4. Given a vertex of r G TM;h,n that is not an endpoint, we can consider the 
subtrees of r with root f , which correspond to the connected components of 
the restriction of r to the vertices w > v. If a subtree with root v contains 
only V and one endpoint on scale + 1, it will be called a trivial subtree. 

5. With each endpoint v we associate a factor V{tp^'^'''''^ + tp^^'^'^^) and a set x„ 
of space-time points (the corresponding integration variables in the x-space 
representation) . 

6. We introduce a field label f to distinguish the field variables appearing in 
the factors V('?/'^*''' -* + ip^^'^^^) associated with the endpoints; the set of field 
labels associated with the endpoint v will be called J^; note that if v is an 
endpoint |Ji,| = 4. Analogously, if v is not an endpoint, we shall call 
the set of field labels associated with the endpoints following the vertex v; 
x(/), e{f), a{f) and p(/) will denote the space-time point, the e index, the 
a index and the p index, respectively, of the Grassmann field variable with 
label /. 

In terms of trees, the effective potential V'-^^ < /i < M (with V^^\i)^'^'^''^) 
identified with Vo(^(*-'^-))), can be written as 

oo 

VW(^(--)+^[i.^l) + /3|A|e;,+i = ^ VW(r,V^('-''-)+V^[^'''l) , (1126) 

where, if Vq is the first vertex of r and Ti, . . . , (s = St,o) are the subtrees of r 
with root fo, V(^)(r, ■i/''^*-'"-) -|- is defined inductively as: 

VW(r,V'[°'^l) = tll!^f,^^,[v(^+i)(ri,^[°''^+i]);...;V('^+^^ . 

(16127) 

where ?/'[°'''] := V^^*-''-) + ^[^''^^ and, if r is a trivial subtree with root on scale M, 
then VW(r,^[0'^^]) = V(^[°'*^l). 

For what follows, it is important to specify the action of the truncated ex- 
pectations on the branches connecting any endpoint v to the closest non-trivial 
vertex v' preceding it. In fact, if r has only one end-point, it is convenient to 
rewrite V^'^H^, ^'°'''') = £l+i£l+2 " ■ ■ Sliiy^il^"^^^^^)) as: 

VW(r, V;[°''^l) = V(^[°''^l) + £l^, ■ • ■^i;(V(^[°'^^]) - V(^[°''^l)) . (1128) 
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Now, the key observation is that, since V(V') is defined as in Eq. fl4[lll) and since 
our exphcit choice of the ultraviolet cutoff makes the tadpoles equal to zero (i.e., 
[g^^\0)]pp = 0), then the second term in the r.h.s. of Eq. (161281) is identically zero: 



^.^+i---^J/(V(V^f°'*^l)-V(^[°''^l))=0 



P29) 



for all < /i < M. Therefore, it is natural to shrink all the branches of r G TM;h,n 
consisting of a subtree r' C r, having root r' on scale h' G [h, M] and only one 
endpoint on scale M + 1, into a trivial subtree, rooted in r' and associated to 
the factor V{il)^^'^'^). By doing so, we end up with an alternative representation 
of the effective potentials, which is based on a slightly modified tree expansion. 
The set of modified trees with n endpoints contributing to V^'*-' will be denoted 
by TM-h,n', every r G Tu-Xn is characterized in the same way as the elements of 
TM-h,n, but for two features: (i) the endpoints of r G Tu-Xn are not necessarily 
on scale M + 1; (ii) every endpoint w of r is attached to a non-trivial vertex on 
scale — 1 and is associated to the factor V{il)^'^'^''~^^). See FiglBl In terms of 




h h-l-l 



M M+1 



FIG. 6. A tree r € TM;h,n with n = 9: the root is on scale h and the endpoints are on 
scales < M + 1. 

these modified trees, we have: 



n=l TeTM:h,n 



(030) 



where 



VW(r,^[o,/.]^ 



\s+l 



S\ 
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and, if r is a trivial subtree with root on scale k G [h,M], then V^'^^t,^^^'''^) = 

Using its inductive definition Eq. fl6[3ip . the right hand side of Eq. fl6[30p can 
be further expanded (it is a sum of several contributions, differing for the choices 
of the field labels contracted under the action of the truncated expectations S'^^ 
associated with the vertices v that are not endpoints), and in order to describe 
the resulting expansion we need some more definitions (allowing us to distinguish 
the fields that are contracted or not "inside the vertex v"). 

We associate with any vertex v of the tree a subset Pt, of 1^, the external fields 
of V. These subsets must satisfy various constraints. First of all, if v is not an 
endpoint and Vi,...,Vs^ are the s^, > 1 vertices immediately following it, then 
Pv ^ UjP^.; if V is an endpoint, = Iv If v is not an endpoint, we shall denote 
by QiJi the intersection of P^ and P^,.; this definition implies that P^ = VJiQ^-. The 
union of the subsets P^,- \ Q^^ is, by definition, the set of the internal fields of 

and is non empty if > 1. Given r G TM-h,n, there are many possible choices 
of the subsets Py, G r, compatible with all the constraints. We shall denote 
by Vt the family of all these choices and by P the elements of Vr- With these 
definitions, we can rewrite V^^\T,-ip^'^'^^) as: 

VW(r,^[°'^l)= 5^ V(^)(r,P), 

vW(r,p) = J rfx,,v^g:yi5^^)(x.o) , mm 

where = U/e/Jx^}, 

f&Pv 

and Kl'^p^\:Ky^^) is defined inductively by the equation, valid for any t> G r that 
is not an endpoint, 

1=1 

where "^Ap "\q has a definition similar to Eq. fl6[33p . Moreover, if Vi is an endpoint 

Kv'^'"~^^\-Xy.) = U; if Vi is not an endpoint, Kv'^'"~^^^ = kI'I^'^^\ where Pj = 
{PuijW G Tj}. Using in the r.h.s. of Eq. ( 16134]) the determinant representation of 
the truncated expectation discussed in the previous section, we finally get: 

TeT TeT 
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where 

W^r,P,T(x.J = P36) 



f/"| n ^ / ^^^T„(t.)detG"-^"(t,) 



V not 
e.p. 



and (^'^"'"^"(t^) is a matrix analogous to the one defined in previous section, with 
g replaced by g^'^^ Note that Wr,p,T and, therefore, V^'^\t,P) do not depend 
on M: therefore, the effective potential V*-'*-' (ip) depends on M only through the 
choice of the scale labels (i.e., the dependence on M is all encoded in TM;h,n)- 
Using Eqs.f l6[35]) -f l6[36p . we finally get the bound: 



-i-|dxi...cix2d<Uxi,...,x20i< if^r E EE- 

^' ' n>max{l,i-l} reTM-hn P&Vr Te"^ 

■ f JJd(x(£) -x'(£)) Y[ ^max|detG'^'"^"(t,)| JJ Wg^^'-^^ii) - 



V not " ieT^ 
e.p. 



(16137) 

Now, an application of the Gram-Hadamard inequality Eq.f l5[6p . combined with 
the representation Eq. (l6ll6p and the dimensional bounds Eq. (l6ll8p . implies that 
|detG'^-^"(t,)| < (const.)S>=i |p.J-|p.|-2{s„-i) _ ^gg^ 

By the decay properties of (^'■'^•'(x) given by Eq. (l6ll5p . it also follows that 

V not ^ feTi, V not '"' 

e.p. e.p. 

(039) 

Plugging Eqs. fl6T38|) - fl6T39|) into Eq.f leTBTj) . we find that the l.h.s. of Eq. (l6T371) can 
be bounded from above by 

E E E E(^°^^t-)"if/r[n 7t2"'"^""'^] • (040) 

n>max{l,Z-l} reTii/.fe „ P&Vt TeT v not '"' 

' ' \Pvq\=21 e.p. 

Using the following relation, which can be easily proved by induction, 

Y hv{sv -l) = h{n-l)+Y (hv - K>){n{v) - 1) , (041) 



V not V not 

e.p. e.p. 



where v' is the vertex immediately preceding v on r and n{v) the number of 
endpoints following v on r, we find that Eq. fl6[40p can be rewritten as 



n>max{l,«-l} reTAf h „ PePr TeT v not 

' ' \Pvo\='2l e.p. 



E EE E(^o^«t-)"it^r2-'^^""'^f n 

(042) 
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where, by construction, we have that n{v) > 1 for any vertex v of t E TM;h,n 
that is not an endpoint (simply because every endpoint f of r is attached to a 
non-trivial vertex on scale — 1, see the discussion after Eg. (16 1 291) ). Now, the 



23 



number of terms in XItgt '^^'^ bounded by (const.)" s„! (see 

Appendix A. 3. 3]); moreover, \P^\ < An(v) and n{v) — 1 > max{l, ^^}, so that 



n{v) — 1 > I + Therefore, 



/3|A| 



J cixi---rfx2z|iyl,'^2jxi,...,X20| < Yl (const.r|f/r2-'^("-') ■ 

?i>max{l,i— 1} 

• 5^ ( n 2-^('^-V)) ( n 2-1^^'!/'') . (043) 

t&Tm h n"" not e.p. PgVt not 

|P„(,|=2Z e.p. 

Now, the sum over P can be bounded as follows: defining py := |P„| (so that 
Pv <Pv:,^ ^PvsJ, then 

P^Vt V not 1) not p„ \ / 

The sum in the r.h.s. can be computed inductively, starting from the root and 
moving towards the endpoints; in fact, at the first step 

Y 2"^"°/'' ^ " ' ^ ^'^"^ ) = + 2-1/1^)^"^ ; (1145) 

at the second step. 



j|(l + 2-Vi6)P., ^2-^"./i6 (^'^^■'^ ^ ■ ■ ■ ^^''-"^ ] = \[ J](l + 2-ii+2 
i=i p„^. \ P-^j ) j=ij'=i 

(16146) 

so that, iterating, 

n E^-"'" h^''^^"'-] = n [E2--r 047) 

V not p„ \ / « e.p. fc=l 

e.p. 

where is the distance from v to the root and we used the fact that p^ = 4 for 
all endpoints v. Plugging Eq. (l6[47p into Eq. (l6[44p . we get 

E ( n < (^171^)" . ei48) 

PeVr ^ not 

Similarly, one can prove that 

( n s-if-M) < , 6,49) 

r€TM;h,n ^ not 
e.p. 
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uniformly in M as M — )■ oo, see 23|, Lemma A. 2]. Collecting all the previous 
bounds, we obtain 



I dxi---rfx2HW^lJi/>i,...,X20| < Yl (const.)1t/r2-'^("-i) , 

ri>max{l,i— 1} 

(Elso) 

which implies Eq.f l6[10]) . The constant Ch can be bounded by the r.h.s. of 
Eq. fl6[37l) with / = and n > 2 (because the contributions to with / = 1 
are zero, by the condition that the tadpoles vanish), which implies 

Ch < ^(const.)"|f/r2-''^"-^) < (const.) If/lV^ . PSI) 

n>2 

Therefore, Fq^m = X^fcii^fc is given by an absolutely convergent power series in 
U, as desired. A critical analysis of the proof shows that all the bounds are 
uniform in M, /3, A and all the expressions involved admit well-defined limits as 
M,/3, |A| —7- oo. In particular, this implies that Fq = limM^oo -fo,M is analytic in 
U (and so is its /3, |A| — )■ oo limit, which is reached uniformly) in the complex 



domain \U\ < Uq, for Uq small enough. See [2J, Appendices C and D] for details 



on these technical aspects. This concludes the proof of Lemma Ell. 



7. THE MULTISCALE INTEGRATION: THE INFRARED REGIME 

We are now left with computing 

Sm,/3,a = e-^l^l^° I P(#(--))e-^«(^"-^-') . dZll) 

We proceed in an iterative fashion, similar to the one described in the previous 
section for the integration of the large values of k^. As a starting point, it is 
convenient to decompose the infrared propagator as: 

g^-^-\^ -y) = J2 e'-^'-^'-'^gl^'K^ - y) , (02) 

UJ=± 

where, if k' = (/cq, k'), 

^i^°Hx-y) = -^ E Xo(|k'|)e--'(^-) f ^I'l -^*Ck'+PF)\ " 

(133) 

and B'^^L = B'^P X Bt, with B^ = {fGi + f G2 - Pf , < ni,n2 < L - 1}. 
Correspondingly, we rewrite ■?/;(*■'■■) as a sum of two independent Grassmann fields: 

Ul = ± 
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and we rewrite Eq. (l6[6p in the form: 



(05) 

where V(°)(V'(-°)) is equal to Vo(V'^'-''-^), once ?/;(*-^-) is rewritten as in fl7[4p . i.e., 
V(°)(^^-°^)= (06) 



a;i,...,aJ2n=± 



5:(/;iAi)-- E E E nc:^. 



n=l 



cri,...,o-„=t4 pl,.--,P2n = l,2 k'^,...,ki,^ j = 



211 



i=i 



j,P2j-l,'^2i-l'^X2j,o-j, 

n=l q_,p,ui j=l 



P2j ,^2j 



with: 

1) 0; = (wi, . . . ,a;2n), £ = (o-i, • • • and p% = (0,^1"); 

2) W^Sp,.(k;, . . . , k'2„_i) = W^M,2n,p(k; + . . . , k'2„„i + p^-^), see (EE 

3) the kernels W^2°|p_(^ (xi , . . . , X2„) are defined as: 



(137) 



2n 



(/3|A|)-- 5^ e^^.-(-)^''^-wit,Jk;,...,k',^^^ . 



k' k' 



Moreover, P^^Q^Aoidil'^-'^^) is defined as 



^Xo,Ao(#^-°^)=Aro-l 



Xo(|k' 



)>0 



n n 



f7=nxo(|k'|)>0 



(138) 



exp { - m\)-' E Xo ^(iki)^g°!.:.Ao,.(k')^i5.;.} 



i^=± k'eB' 



/3,-L 



where: 



-t;ofi*(A;'+pJ^) 



-<o^o + so(k') fo(i/i;i - ujk'^) + to,a;(k') 
^t;o(-2A;'i - uk'^) + ^^^^(k') -^CoA^o + So(k') 

Ao is chosen in such a way that / P^^^Aoid^p^-'^'^) = 1, Co = 1, = and 
|to,.(k')|<C|kf. 
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It is apparent that the ip^-^^ field has zero mass (i.e., its propagator decays 
polynomially at large distances in x-space). Therefore, its integration requires 
an infrared multiscale analysis. As in the analysis of the ultraviolet problem, 
we define a sequence of geometrically decreasing momentum scales 2^*, with h = 
0, —1, —2, . . . Correspondingly we introduce compact support functions /h(k') = 
XoC^'^l^^'l) - Xo(2"''+^|k'|) and we rewrite 



Xo(lk'l) = Yl ^(k') 



(09) 



h=ha 



with hp := [log2 (||)J (the reason why the sum in Eq. fl7[9p runs over h > h 



is that, if k' G B'^^j^ and h < /i^, then //i(k'), simply because |k'| > vr//?). The 
purpose is to perform the integration of (171 5 p in an iterative way. We step by 
step decompose the propagator into a sum of two propagators, the first supported 
on momenta ^ 2^, h < 0, the second supported on momenta smaller than 2^. 
Correspondingly we rewrite the Grassmann field as a sum of two independent 
fields: ijj^-'^^ = ip^'^^ + ^(-''~^) and we integrate the field In this way we 

inductively prove that, for any h <0, Eq. fl7l5p can be rewritten as 



-/3|A|Fh 



,(# 



(<'i))g-v(''HV''-''') 



dlllO) 



where Fh,Ah,V^^'' will be defined recursively, x/i(|k'|) = Ylk=h /fc(k') and 
^Xh,Ahi'^'^^~'^^) is defined in the same way as Py^f^^Aoid^^-^^) with tp^-^Kxo^ 
Ao,LoXo, Vo,So,to^^ replaced by ^^-^\xh,Ah,uj,Ch,Vh,Sh,th,u^, respectively. More- 
over VW(0) = and 



2n 



n=l 



(011) 



E 

n>l 



i<h)- 

X2j,o-j,P2j,ti^2j 



X2ri 



Note that the field 'ip^'^^, whose propagator is given by x/i(|k'|)[y4if'^(k')]^\ 
has the same support as Xh, that is on a neighborood of size 2^ around the 
singularity k' = (that, in the original variables, corresponds to the Dirac point 
k = Pf)- It is important for the following to think 1^2^,p,a;5 ^ < 0, as functions 
of the variables {Ck,Vk}h<k<o- The iterative construction below will inductively 
imply that the dependence on these variables is well defined. The iteration con- 
tinues up to the scale hp and the result of the last iteration is 'Bm^p^a = e~ 
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Localization and renormalization. In order to inductively prove Eq.f l7[10l) we 
write 

where 

, uj=± Xh{|k'|)>0 

and 7^V('*) is given by Eg. (171 lip with replaced by X]^2' ^'^^^ contains 

only the monomials with four fermionic fields or more. The symmetries of the 
action, which are described in Appendix [C] and are preserved by the iterative 
integration procedure, imply that, in the zero temperature and thermodynamic 
limit, 

^^'i''''^'' 

for suitable real constants Zh, Sh- Eq. fl7[14p is the analogue of Eq. fl6[12p : its proof 
is completely analogous to the proof of Lemma El 2 and will not be belabored 
here. 

Once that the above definitions are given, we can describe our iterative inte- 
gration procedure for h < 0. We start from Eg. (171100 and we rewrite it as 

^x.,A.(#^-''^) e-^^*'^(^*"'")-^^""(^*"'")-'3|^l^'' , (17115) 
with 

Xh(|k'|)>0 

£VW(^(<^)) = (/3|A|)-i^ J2 ■ (1316) 



k' 



■^k',.,. ys,{-zk[ - uk2) + T*Jk') -zz^ko + a,(k') 1 "^"^''^'^ • 
Then we include CV^'^^ in the fermionic integration, so obtaining 

^x.,A.-i(^^^-'^) e-^^^'^('^^"")-^l^l(^'^+^'') , (017) 
where 

= ^EEE^Tr{[x.(kOA-(ko<';L,.)(k')]"} mis) 

^' ' U],a k' n>l 
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is a constant taking into account the change in the normahzation factor of the 
measure and 



(dig) 

with: 



C^_i(k') = a + ZhXh{\^) , ^^h-i(kO =Vh + 6hXh{^^) , (020) 

s,_i(kO = s,(kO + a^(k')xh(k') , t„-i,.(k') = t,,^(k') + r,,^(k')xh(k') . 

Now we can perform the integration of the il)^'^'^ field. We rewrite the Grassmann 
field if)'^-^'^ as a sum of two independent Grassmann fields ■^(-'^~^) + ■?/)('^) and 
correspondingly we rewrite Eq. fl7[17p as 



(1321) 



where 

Aft-l,a.(k') 



-iCh-iko + s/j_i(k') Vh-i{ik[ - uk'^) + t/i_i,a;(k') 
t;/,_i(-iA;; - w/c^) + tLi,a;(k') -<fc-iA;o + Sfe_i(k') 

(17122) 

with: 



Ch-i = Cfc + 2;/! , v/i-i =Vh + Sh , 

Sh-i{k') = s,,(k') + ah{k') , t;.-i,.(kO = U,^(k') + r,„^(k') . ((7123) 

The single scale propagator is 

/ = - x^) , (1324) 



where 



eHxi-X2) = -^ 5^ e-^'^'(---)A(k')k_i,.(k')l"'. (1325) 



/3|A| 



After the integration of the field on scale h we are left with an integral involving 
the fields ■?/)(-'*~^) and the new effective interaction V^'*"^-', defined as 



with V(''"^)(0) = 0. It is easy to see that V^^"^) is of the form Eg- flTTTT]) and that 
Fh-i = Fh + eh + eh- It is sufficient to use the identity 

n>l 
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where £'^{X{ip^'^^)', n) is the truncated expectation of order n w.r.t. the propaga- 
tor which is the analogue of Eq. fl6[7p with ^/)("-*'-) replaced by i})'^^'' and with 
Pi^dip'^u.v.)^ replaced by P^^^^^_^(#('^)). 

Note that the above procedure allows us to write the effective constants {(h, Vh), 
h < 0, in terms of iCk,Vk), h < k < 0, namely Ch~i = (3i{{Ch,Vh), • • • , (Co,^^o)) 
and Vh-i = l3l(y{C,h^Vh), • • • , (Co?'^^o))) where /3* is the so-called Beta function. 

An iterative implementation of Eg. (17127]) leads to a representation of V^'*^ 
/i < 0, in terms of a new tree expansion. The set of trees of order n contributing 
to V^^''\ip^-'^'>) is denoted by 7h,n- The trees in Th^n are defined in a way very 
similar to those in TM-h,n, but for the following differences: (i) all the endpoints 
are on scale 1; (ii) the scale labels of the vertices of r G Th,n are between h + 1 and 
1; (iii) with each endpoint v we associate one of the monomials with four or more 
Grassmann fields contributing to TZV^'^\ip^-^^) , corresponding to the terms with 
n > 2 in the r.h.s. of Eq.( ]7l6p (here TZ is the linear operator acting as the identity 
on Grassmann monomials of order 4 or more, and acting as on Grassmann 
monomials of order or 2). In terms of these trees, the effective potential V^'*^ 
h < —1, can be written as 

oo 

VW(V;(^'^)) + /3|A|efc+i = J2J2 V(^)(r,^(^^)) , (028) 

where, if Vq is the first vertex of r and ri, . . . , (s = s^q) are the subtrees of r 
with root Vq, 

VW(r,V^(2^'')) = ^zll!li^T^^[7^v(^+l)(rl,^(^^+l)); . . . ; 7^V('^+l)(r„ , 

(17129) 

where, if is trivial, then /i + 1 = and 7^V(°)(r„ V^^^'^+i)) = 7^V(°)(^(^''+^)). 
Repeating step by step the discussion leading to Eqs. fl6[32p . (l6[35p and f l6[36p . and 
using analogous definitions, we find that 

with 



.i=l 



\\{—\ I ^^T„(t.)detG^"'^''(t„) 

I V not J 
e.p. 



'a;W,<^'W5<xW,a'W [c/2;}(x(£) - X (£))] 



(1331) 
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and: v*, i = 1, . . . , n, are the endpoints of r; kII^'\-x^*) is the kernel of one of the 
monomials contributing to TZV^^''{tp^-'^''^); G^'^ is a matrix with elements 

^/J = hf)Af')^'^{f)Mf')^'^U)Mf') - ^(•^'))]p(/),p(/') ' (022) 

Once again, it is important to note that G^''^ is a Gram matrix, i.e., defining 
e+ = e-|- = (1, 0) and e_ = = (0, 1), the matrix elements in Eq. fl7[32p . using a 
notation similar to Eq.f l5[7p . can be written in terms of scalar products: 

GjJ, = (17133) 

= {ui(j) (g) e^(/) (g) e^(/) (g) Ax(/),p(/),<^(/), Ui(//) (g) e^(/,) (g e^(/,) (g Bx(/'),p(/'),a;(/')) ) 

where 



^ "'"^ /3|A|,,|:iJdetZ._i,.(k')|V4 ' 



Pi* ' 



1 „-jk'(z-x) /TTiTTT 

Kh-ih - s/i-i(k') Vh-i{ik[ - uk'2) + U-i,<^(k') 

Uft_i(-iA;; - w^a) + iLi,^(k') </x-i^o - s/,_i(k') 

so that, using that Sh-i is purely imaginary, 

for a suitable constant C. 

Using the representation Eqs. fl7[ 28 p . ( I7[ 29 p . (171 3 1 p and proceeding as in the 
proof of Lemma Ol, we can get a bound on the kernels of the effective potentials, 
which is the key ingredient for the proof of Theorem 12.11 and is summarized in 
the following theorem. 

Theorem 7.1 There exists a constants Uq > 0, independent of M, /3 and A, such 
that the kernels W^2i'|paj(^i' • • • 5^2/) in Eq. ^7\ll\ ), h < —1, are analytic functions 
ofU in the complex domain \U\ < Uq, satisfying, for any Q <6 <1 and a suitable 
constant > (independent of M, /3, A), the following estimates: 

^ / rfxi ■ ■ ■ d^2i\W^lj^,, . . . , < 2'^'-''^'^ (Co |f/|)— (I'^-i) . (17136) 

Moreover, the constants Ch andch defined by Eq.^J~W^ and Eq.^.27) are analytic 
functions of U in the same domain \U\ < Uq, uniformly bounded as |e/j| + |e/j| < 
C0\U\2^''^^^\ Both the kernels W^^fp^l^i, . . . ,X2;) and the constants ChjCh admit 
well defined limits as M, /3, |A| — > 00, which are reached uniformly (and, with 
some abuse of notation, will be denoted by the same symbols). 
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Before we present its proof, let us show how Theorem 17.11 imphes, as a corollary, 
Theorem 12.11 It is enough to observe that, by Proposition 14.11 and by the multi- 
scale integration procedure described above, fpiU) = Fq + Y^h=hfM^ + ^h), with 
-fo, Gh^e-h analytic functions of U for U small enough (see the discussion at the end 
of Section |6] and the statement of Theorem 17. ip . Since \eh\ + \eh\ < Ce\U\2'^'^'^^^\ 
the sum Y^h=hpi^h + e-h) is absolutely convergent, uniformly in hp] therefore, both 
f/3{U) and its /3 — oo limit are analytic functions of U for U small enough. This 
concludes the proof of Theorem 12. 1[ 

Proof of Theorem 17. IL Let us preliminarily assume that, for h' < h < —1, 
and for suitable constants c, c„, the corrections Zh, Sh, o"/i(k') and T/i(k') defined in 
Eq. (l7[14l) and Eq. (l7[16p . satisfy the following estimates: 

max {\zhl\Sh\} < c\U\2'^ , (037) 
sup UKMrnAK'ThA^'m < c„|f/|22('^'-'^)2(i+^-")'^ . 

2h'-i<|k'|<2'»'+i 

Using Eg. (171371) we inductively see that the running coupling functions (h, Vh, 
s/i(k') and th(k') satisfy similar estimates: 

max {\a - l\, \vh- Vol} <c\U\ , (038) 
sup {PJ^,s,(k')||, ||9^,(t/.,.(k') - to,.(k'))||} < c„|f/|22('^'-^)2(i+^-)'' . 

2h'-i<|k'|<2'i'+l 

Now, using the definition of g!;^\ see Eq. (l7[25p and Eq. (l7[19p . and the bounds 
Eq. (l7[38p . we get (proceeding as in the proof of Eq. flAISp in Appendix Rl). 

92/1 

for all i^' > and a suitable constant Ck- Using the tree expansion described 
above and, in particular, Eqs. (171 28|) . (I7[30p . (1713 ip . we find that the l.h.s. of 
Eq. (l7[36p can be bounded from above by 



n 



.1=1 



(17140) 



\Pvo\='2l 



n ;^max|detG'^-^"(t.)| J] Wal'iW)-^'^ 
V not " £er„ 

c.p. 

where T* is a tree graph obtained from T = U^T^, by adding in a suitable 
(obvious) way, for each endpoint v*, i = 1, . . . ,n, one or more lines connecting 
the space-time points belonging to x„f . 
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An application of the Gram-Hadamard inequality Eq. fl5[6p . combined with the 
representation Eq. (171331) and the dimensional bound Eq.f l7l35p . implies that 

|detG"^-^''(t,)| < (const.)^-il^"'l-l^''l-'(^"-^)- 2'^"(^£i l^-^l-l^^l-^^^^-^)) . I^Al) 

By the decay properties of g^\-x), Eq. (171 39 p . it also follows that 

n ^ / n ^(^w - ^'w) - < n ^2-^^(--i) . 

V not ^' l^Ty V not ^' 

e.p. e.p. 

(17142) 

The bound Eq. (l6[10p on the kernels produced by the ultraviolet integration im- 
plies that 

/n n 
n «£) -x'(£))nii^i^Hx<)i < wc^m'i-' , (17143) 

teT*\U„T„ i=l i=l 

where pi = \Pv*\- Combining the previous bounds, we find that Eq. (l7[40p can be 
bounded from above by 

n>lTeThn PePr TeT not i=l 

' \Pvo\=2l e.p. 

(17144) 

Let us recall that n{v) = J^i-v'yv ^ number of endpoints following v on r, 

that v' is the vertex immediately preceding f on r and that is the number of 
field labels associated to the endpoints following t> on r (note that |/,;| > 4n(f)). 
Using the following relations, which can be easily proved by induction, 

V not e.p. i=l v not e.p. 

K{s, - 1) = h{n - 1) + {K - K>){s, - 1) , (17145) 

V not e.p. V not e.p. 

we find that Eq. (l7[44p can be bounded above by 

Y ^ C'"2^(^~l^"ol+l^''ol-3'') ■ 

\P.o\=2l 

1 " 



V not i=l 
e.p. 



Using the following identities 



2hn 2('^^-'^v'M^) = 2^"' , ([7147) 



V not V e.p. 

e.p. 



2^11-^1.01 JJ^ 2^^""'^^'^^^^^ = Y[ 2'^"'''^"' , ([7l4J 



V not V e.p. 

e.p. 



41 



combined with remark that all the endpoints are on scale 1 (so that the right 
hand sides of Eqs. (171471) - f l7l48p are equal to 1), we obtain 

1 



(3\A\ 



n>i T&Th n PePr TeT 



1 

1=1 



(17149) 



not 
e.p. 



Note that, if v is not an endpoint, |P^| — 3 > 1 by the definition of TZ. Now, note 
that the number of terms in XItst "^^^ bounded by C"' Hi; not ep -^f'- Using 
also that |P„| - 3 > 1 and |P^| - 3 > |P„|/4, we find that the l.h.s. of Eq. fTMl 
can be bounded as 
1 

W\ 



i-/rfx,...rfx,i<)jx,,...,x,oi<2'^(^-i"^^^ y: ■ 

n 



(II150) 



i; not 
e.p. 



FgVt ^ not 



i=l 



Proceeding as in the previous section, we get the analogue of Eq. (l6[48p : 



PGPt f not 
\Pvo\=-2l e.p. 



i=l 



Finally, using that 2-^^^-"-^^') < 2^^ and that, for < < 1 (in analogy 

with Eq.dnBnD), 

T&Th,n not e.p. 

and collecting all the previous bounds, we obtain 

^ /c^xi ■ ■ ■rfx2d<^,^(xi, . . . ,X2z)| < 2^^(3-1^-1+^) 5^ C"1[/r , (051) 



/3|A| 



n>l 



which is the desired result. 

We now need to prove the assumption Eq. (l7[37p . We proceed by induction. It 
is easy to see that the assumption is valid for h = 0; in fact, 

.(0) 
2 

(052) 



|;^0| = |5fco<d,l),(.,.)(PF)l , l^ol = |5,,#2,(1,2).(.,.)(P^^)| 

\ao{k')\ = \y2 ds'k'^Kd^d^W2,{i,i),(uj,uj){s'k' 



„I/=02 



^0 



ro(k')| = \Y^ / ds ds'k'^Kd,d,W2Xi,2U^,u.){s''k:) 
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so that, 



max{|2;o|, |5o|} < y rfxc/y ||x - y|| ||iy2,(i^,L^)(x - y)|| , (053) 

sup {||a^,ao(k')||,||ai^,ro(k')||} < -rr-r / rfx^;y||x-y|r2||1^2,(.,.)(x-y)|| , 

which can be bounded via the same strategy as the proof of Eg. (161 101) . yielding 
Eg. (1 71371) for /i = 0. Similarly, assuming that Eg. fl7[37p is valid for all /i > /c + 1, 
the guantities of interest ioi h = k can be bounded by 

max{|z,|,|4|} < ^/ rfxdy||x-y|| ||<(L,.)(x, y)|| , (054) 



and 



sup {PJ:,a,(k')||,PJ:,r,,^(k')||}< 

-i<|k'|<2'''+i 



C2 



2h' 



The same proof leading to Eg. (1 715 II) shows that the r.h.s. of Eg. fl7[54p can be 
bounded by the r.h.s. of Eg.f l7[51]) times 2~^ (that is the dimensional estimate 
for ||x — y||), and that the r.h.s. of Eg.f l7[54|) can be bounded by the r.h.s. of 
Eg. fl 715 II) times (const. )"2^'^'2~'^"+^)'^ (where 2"'^("+^) is the dimensional estimate 
for ||x — y||"^^). 

It remains to prove the estimates on eh,eh- The bound on eh is an imme- 
diate corollary of the discussion above, simply because can be bounded by 
Eg. (171 40 p with / = 0. Finally, remember that Ch is given by Eg. (l7ll8p : an explicit 
computation of A^l{k')W^'^^^^^{k') and the use of Egs. (171371) -( 171381) imply that 

Pi(k')<(i,.)(k')ll < C\U\2<^\ from which: |e.| < C'2^^j:n>iiC\U\2<^'r, as 
desired. This concludes the proof of Theorem 17.11 and, therefore, as discussed 
after the statement of Theorem 17.11 it also concludes the proof of analyticity of 
the specific free energy and ground state energy. ■ 



8. CONCLUSIONS 

In conclusion, I presented a self-contained proof of the analyticity of the specific 
free energy and ground state energy of the 2D Hubbard model on the honeycomb 
lattice, at half-filling and weak coupling. The proof is based on rigorous fermionic 
RG methods and can be extended to the construction of the interacting correla- 
;ions, i.e., the off-diagonal elements of the reduced density matrices of the system 



2J], and to the computation of the universal optical conductivity ^8|. Such con- 



struction shows that the interacting correlations decay to zero at infinity with 
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the same decay exponents as those of the non-interacting case. The "only" effect 
of the interactions is to change by a finite amount the quasi-particle weight Z^^ 
at the Fermi surface and the Fermi velocity v. 

The example presented here is the only known example of a realistic 2D in- 
teracting Fermi system for which the ground state (including the correlations) 
can be constructed. The main difference with respect to other more standard 
2D Fermi systems is the fact that here, at half-filling, the Fermi surface reduces 
to a set of two isolated points. This fact dramatically improves the infrared 
scaling properties of the theory: the four-fermions interaction, rather than being 
marginal, as in many other similar cases, is irrelevant; this is the technical point 
that makes the construction of the ground state possible and "relatively easy" . 

It is natural to ask how the system behaves in the presence of electromagnetic 
interactions among the electrons, which is the case of interest for applications to 
clean graphene samples. In this case, the system has many analogies with (2 + 1)- 
dimensional QED. The four-fermions interaction, rather than being irrelevant, is 
marginal, and the fixed point of the theory is expected to be non-trivial. The long 
distance decay of correlations is expected to be described in terms of anomalous 
critical exponents and the effective Fermi velocity is expected to grow up to the 
maximal possible value, i.e., the speed of light. The specific values of the critical 
exponents suggest that local distortions of the lattice (in the form of the so- 
called Kekule pattern ^o|) are amplified by electromagnetic interactions: this 
led us to propose a possible mechanism for the spontaneous formation of the 
Kekule pattern, via a mechanism analogous to the ID Peierls' instability 27 1. 
All these claims have been proved so far only order by order in renormalized 



perturbation theory 26|, |27(|; proving them in a non-perturbative fashion is an 



important open problem, whose solution would represent a corner stone in the 
mathematical theory of quantum Coulomb systems. 



Appendix A: Dimensional estimates of the propagators 

In this Appendix we prove the dimensional bounds on the propagators used 
in Sections miMEHT] and, in particular, the bounds \\g\\oa < CM and \\g\\\ < C(3, 
stated right before Eq.( l4[26p . and the estimates Eqs.f l6[15|) .f l7[39p . The basic idea 
is to decompose the propagator in Eq. fl4[5p . 

'^' KZ^-n ki + vS\Sl{k)\^ \-v„a(k) tk„ I 
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as a sum of single scale propagators: 

M 

^(x) = Y.g^^\^) + ^i'Hx) , (02) 

h=l a;=ib h=hp 

where: (i) if 1 < < M, then g'^'^^x.) is defined as in Eq. fl6[14p : (ii) if /i^ : = 
Llog2 (if )J and < < 0, then 

with fh(k') = xo(2"^|k'|) - xoi2'^^%'\)- The main issue is to prove Eq. llGTTHll . 
and the analogue of Eg. (171391) . i.e., 



Ck2 



2h 



where ||x||2 = |xo|^ + |x|i, for all K >0. Note that ||^||oo < CM and ||^||i < C(3 
are immediate corollaries of Eqs.( ]A[4|) -( ]A[5l) . In fact, plugging these bounds into 
Eq.dffl), we find: 

M 
h=l uj=± h=hp 



< C4 V / rfx r^r^ + 2C4 V /" rfx- 



< 



1 + (2^^10:01/3 + \x\a)^ Va) 1 + (2'^l|x||)^ 

AI 

which are the desired bounds. 

So, let us start by proving Eq. (]A[4|) . We denote by 9k the discrete derivative, 
and by 9^ = ■ 9k, with /i G {0, 1,2}, the discrete derivative in the direction 

(here Cq = (1,0,0), ei = (0,1,^) and 62 = (0,|,— ^)), defined as follows: 



given a compact support function -F(k) with k G -BI"^'', we let 



9„nk) = ev.a.F(k) = ^(;i±i4M^. a, 

with AA;o = y and AA;i = Afcg = |f . Note that, if 



F(x) = e-^^-''F(k) , (jAl9) 
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then, defining = x ■ e^, 



keB 



(oo) 
'13, L 



dAlio) 



so that, using the fact that |xo|/3 < f'i/3(a;o) and |x|a < + d1{x2)y^'^ 

(here d/3(a;o) = and diixi) = a^Jgf ) , we find 



l^o|||F(x)|<-4(xo)|F(x)| = - ' 



4 



4 



Ako 



keB 



(oo) 



(011) 



keB 



(oo) 

;3,-L 



keB' 



(oo) 
0,L 



and, similarly, 

l^li|i^(x)|<y(4(5i) + 4fe))|F(x) 



2 ^ 



2 2 ^ gjAfc^Xp — 1 2 



E e-^'^-^F(k) 



keB 



(oo) 
/3,L 



(012) 



y 



El E ^''''"'9lF{k) 



^•=1 keB<-°°J 



< 



37r2 



E K^'i + ^'.mk)!- 



/9,i 



keB 



(<x=) 



Coming back to Eg. flA[4p . recalling that 



Therefore, 



- T 



-ikx 



keB 



(M) 
/3,L 



fu.vXk)Hhiko) ( iko -von*{k) 
kl + vl\n{k)\^ \-von{k) iko 



keB 



(M) 
/3,i 



[A:2 + .;g|fi(fc)|2] 



/i > 2 . 
(013) 

(IA114) 



where in the last inequality we used the fact that, on the support of /„.„.(k)if/i(A;o), 
feg+f Q|f2(/i;) P ~ 2"^^ (here ~ means that the ratio of the two sides is bounded above 
and below by universal constants) and that, moreover, the measure of the support 
of fu.MHhiko) is itself of order 2^ that is (/3|A|)-^ Ek A.-(k)i^h(fco) ~ 2^ 
Moreover, using Eqs. (IA[lip - (IA[12p . we have that for all > and a suitable 
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constant C, 



2N 
ko 



fu.v.(M)Hh{ko) / iko -von*{k) 



/3|A| ^11 La;2 + ^2|^(^)|2 \-t;ofi(A;) zA;o 



fu.v.ik)Hh{ko) ( iko -von*{k) 



^1^1 1^.)" '^ki+vimw\-^om ^ko 



Now, in the first line, every derivative with respect to /cq can act either on the 
denominator k^ + VQ\VL{k)\'^, or on fu.^ (k)Hh{ko), or on the diagonal elements of 
the matrix, iko] in all these cases, recalling that, on the support of fu.v.Hh, \ko\ ^ 
2^, every derivative dk^ can be estimated dimensionally by a factor proportional 
to 2"^^. This leads to the bound 

|a;o|J'^||^7('^)(x)|| < (const. )'^2'' 2-2^^^ 2-^^ , (EllS) 

which differs from Eq.f lA[14l) precisely by the factor 2~^^'^. Similarly, in the 
second line, the derivative with respect to k can act either on the denominator 
fcQ + fQ|f2(fc)p, or on /u.„.(k)if/i(fco), or on the off-diagonal elements of the matrix, 
—vo^{k) or —VQQ*{k); in the first case, dj^ can be estimated dimensionally by a 
factor proportional to 2"'^ while, in the other cases, by an order 1 factor. This 
leads to the bound 

\x\f\\g^''\^)\ \ < (const.)^2'^2-'^ , (016) 

which, if combined with Eqs . (I A [ 1"4|) - (I A [ 1 5"|l . finally implies Eg. (lA[4p . 
The proof of Eg. flA[4p is very similar. In fact, for all > 0, 

II l|27V||-W/ x|| < \\r.2N\ fhiko,k-p^) f iko -Von*(k) 

^ ^" - /5|A| ^f^, i + viMkW V-om ^ko 

(017) 

Now, using the fact that, on the support of fh{ko, k—pp), \ko\, |f^(fc)| ~ 2 , 

we see that every derivative 9k, when acting either on f^, or on the denominator 
ko + "^01^(^)1^5 elements of the matrix, can be estimated dimensionally 

by a factor proportional to 2~^; moreover, the measure of the support of fh is 
itself of order 2^^", that is il3\A\y^J2kfhiko,k-f'p) ~ 2^^. This leads to the 
bound 

||x||2^||^(f)(x)|| < (const.)^23'^2-2^'^2-^ , (018) 



which finally implies Eg. (IA[4p . 
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Appendix B: Truncated expectations and determinants 

In this Appendix we prove Eq. fl5[3p . following j23j, Appendix A. 3. 2]. Given s 
set of indices Pi, . . . ,Ps, consider the quantity £'^{ipp^, . . . , ippj. Define 

Pt = {fe P, : e{f) = ±} Pl) 

and set / = for / e Pf, with i = l,...,|Pf|. Note that E?=i l^/l = 

'^'j=i otherwise the considered truncated expectation is vanishing. Define 

s 

j,j'=l i=l i'=l 

where thf. := ib^, ■ n .^ /. n and, if = , IPJ^I = , |PrL then G is the 
n X n matrix with entries 

SO that 



S (fl'^P^ =detG = jvij exp [- . 



P5) 



Setting X := {1, . . . , s} and 



ip-ii^;i 



i=l i'=l 



we write 



where 



i,i'6X j<j' 



In terms of these definitions, Eq. llBISP can be rewritten as 



£ (^n^^^. j = detG = j Vije-^^^^ . ^9) 

We now want to express the last expression in terms of the functions Wx, defined 
as follows: 

r 

Wx{X,,...,Xr;h,...,tr)= Yl lltk{i)Ve, PlO) 

e£L{X) k=l 

where: 
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1. Xk are subsets of X with \Xi,\ = k and such that 

Xi = {1} , 
-^fc+i 3 Xk ; 



2. L{X) is the set of unordered pairs in X, i.e., the set of pairs {j,j') with 
G X, such that (j, j') is identified with (j', j) and, possibly, j = j'; we 
shall say that £ = is non tivial if j ^ j', and trivial otherwise; 



3. the functions tk{i) are defined as follows: 

tk , if ^ ~ , 



tk{i) 



1 , otherwise 



Pl2) 



where i ~ means that i = {],]') "intersects the boundary" of Xfc, i.e., 
it connects j G with j' ^ X^, or viceversa. See Fig. [3 




FIG. 7. Graphical representation of the sets X^, = 1,2,3. In the example Xi = {!}, 
X2 = {1,2}, X3 = {1,2,3} and ^ = (1,3). The line i intersects both the boundary of 
Xi and of X2, i.e., I ~ dXi and £ ~ 8X2- 



Let us show how to re-express e in terms of the IFx's. The basic step is 
the following: using the definition Eq. flBllOp . we rewrite 



Wx{X^-, ti) = hV{X) + (1 - ti) [V(Xi) + V(X\Xi)] 



Pl3) 



that is, we recognize that Wx{Xi\ ti) interpolates between the full V{X) and two 
of its "proper subsets", V^(Xi), V"(X\Xi). In this way. 



-V(X) 



^-W^{Xrfi)^ / dti 



d_ 



P14) 
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Let us now iterate this construction: let us consider one of the terms in the 
summation in the r.h.s. of Eq.( ]B[14l) : if ^l = (1, j*), we let X2 := Xi U {j*} and 
we note that, by definition, 

Wx{X,,X2;h,h) = t2Wx{X,;t,) + {1 -t2)[Wx,{X,;t,) + V{X\X2)] , Pl5) 

so that 



^-Wx{Xv,ti) ^ ^-Wx(Xi,X2;ti,0) 



-Wx{Xi,X2;h,0) 



d_ 



-Wx{Xi,X2-MM) 



Pl6) 



£2-9x2 

Substituting Eq. dBTTGj) into Eq.f lBTM]) we get: 



Wx(Xi,X2;h,t2) 



g-V(X) ^ ^-WxiX 



^-Wx{XuX2;U,0) _^ 



e.^.i^v. e^^i^v^-^O JO 



Pi7) 



'Wx{Xi,X2;ti,t2) 



ii^dXi £2^8X2 

Using the fact that 

Wx{Xi, . . . , Xp^i, ti, . . . , tp+i) = tp^iWx{Xi, . . . , Xp] ti, . . . ,tp) + 

(1 - Vi) [Wx^^AXi, ...,Xp;t,,...,tp) + V{X \ Xp+O] Pi^ 

and iterating we find 



-ViX) 



nti(4+i) . . .tfc(4+i) e-^-(^-'^^+^^*^'-'*-°) , Pl9) 

where, if r = 0, the summand should be interpreted as equal to e~^^^^^'°\ 
Moreover, by the definition of the WOf 's, if r > 1, 

iyx(Xi, . . . , X,; ti, . . . , 0) = Wx^ (Xi, . . . , X,_i; ti, . . . , + l^(X \ X,,) . 

P2O) 

Using this remark and the notion of anchored tree, which will be defined in a 
moment, Eq. (]B[19P can be rewritten in a more suggestive and convenient way. 
Let T(X) be the set of tree graphs on X, i.e., the set of (|X| — l)-ples of lines 
in L{X) connecting (in a minimal way) all the elements of X. Given a sequence 
of subsets Xi C ■ ■ ■ C Xr as above, we shall say that T G T{Xr) is an anchored 
tree on (Xi, . . . , X^) if its lines can be ordered in such a way that ii ~ (?Xi, 4 ~ 
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(9X2, . . . , ir-i ~ dXr-i- Moreover, given a sequence Xi C ■ ■ ■ C as above and 
a non-trivial line i G -^^(X), we let 

n{e) = max{k : i ~ dXk} , = min{A; : £ ~ dXk} ; ^21) 

if £ is trivial, we let n{i) = n'{i) = 0. Using these definitions, we rewrite Eq. flB[19p 
as: 

r=l XrCX X2,.-.,^r-l T on (Xi,...,Xr) teT 

rr-l 







/"'dti.. 




/o 





g-W^X,(Xl,.--,-^r-i;tl,.--,tr-l) g-y(X\X,) 



where "T on (Xi, . . . , X^)" means that T is an anchored tree on (Xi, . . . , X^). 
Defining 

X2,...,Xr-i T on {Xi,...,Xr) ieT 



,Xr-\;ti,...,tr-\) 



Eq.f lB[22|) can be further rewritten as 

^-v(x) ^ K{Y) e-^(^\^) , P24) 



YCX 
YB{1} 



and, iterating. 



m=l (yi,...,y™) i=l 

where the second summation runs over partitions of X of multiplicity m, i.e., over 
m-ples of disjoint sets Y^, . . . ,Y"^ such that W^^Y"^ = X. Plugging Eq. (lB[25p 
into Eq. (lB[9p gives 

n^^. = E E (-l)^^'"(-ir / Vi;y.l[K{Y^) , P26) 



(/),p{/) 



and f-l)'" 



where V^y^ = H.-^y. [n/eP+ d^x(/),.(/),p{/)^ ^ 
is the sign of the permutation leading from the ordering of the fields in Dip to 
the ones in Wi'Dipy^- In Eq. ( 1B[26"1) . each factor K{Y^), after small manipulations 
of its definition, Eq. (IB[23|) . can be rewritten as 

K{y') = E E \{v,fdti...C dt\y.\_i ■ 

^ ' 2' ' \Y'\-1 

compatible with T 

■ Ylitn'iE) ■ ■ ■t„(,)_i)e-^^"(^')*"'w-*"(^)^^ , P27) 
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where: (i) Y-^ := {mm{j : j G Y"^}}; (ii) the second summations is over sequences 
of subsets Y^ of Y' with \ Y^\ = q and such that Y^ C Y^ C ■ ■ ■ C Y^\.,^_^ C Y^y,^ = 
and T is anchored on (Y"/, . . . , Y|yi|); (iii) in the second hne, if n{i) = n'{£), 
the factor tn'(i) ■ ■ ■ tn(£)~i should be interpreted as equal to 1; (iv) in the exponent, 
if £ is trivial (and, therefore, n\t) = n{i) = 0), to should be interpreted as equal 
to 1. Now, using the analogue of Eq. fl4[20p in a case, like the one at hand, where 
the monomials ip do not necessarily all commute among each other, we have that 
(denoting the elements of Y^ as Y^ = . . . , j|y.|}) 

^ ( n ) = E E i-'^r^'^i^p,. , . . . , J ■ ■ ■ s^i^p,^ , . . . , v^p,^^^ ) 

\i=i / m=i (yi,...,y™) ^ l'^^ 

P28) 

where (— 1)*^ is the parity of the permutation leading to the ordering on the r.h.s. 
from the one on the l.h.s.; note that a is the same as in eqn( ]B[26l) . Comparing 
Eq. flBl28l) with Eq. flBl26l) and using Eq. flBl27l) . we get: 

£^(^P„...,^pJ = (-1)^+^ Yl l^i^U^^I rf^T(t)e-^(*) , P29) 



Ter(x) " eeT 



where we defined: 



s-l 



dPT{t)= Y[{tn'ii)...tn{i)-l)Y[dtg (P30) 



X2,...,Xs-i ieT q=l 
compatible with T 



and 



\/(t)= J2 tn'(l)---tn[e)V,. P31) 

ieL{X) 

Finally, if we use the definitions Eqs. (1B[6l) . (IB[8l) and explicitly integrate the Grass- 
mann variables along the lines of the anchored tree, we end up with 

(^P„ . . . , V^pJ = "^n%./l / / d^T(t)e-^*(*) , P32) 

where V = {Pi, . . . , Pg} and T{V) is the set of anchored trees between the "boxes" 
Pi, . . . ,Ps, i.e., the graphs that become trees if one identifies all the points in the 
same "clusters" Pi (note that now the lines of an anchored tree T G T{V) are 
pairs of field variables {f,f'), rather than pairs of indices {],]'))■ Moreover, ap 
is a sign (irrelevant to the purpose of the bounds performed in this paper), f}, 
f1 are the two field labels associated to the two (entering and exiting) half-lines 
contracted into £, and 

V*{di^) = n ^€(/),.(/),p(/) , = E ^-'W • • • (033) 

/eP e€L(x) 
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where 



VL = { ^ ' _r ^ ^' ' P34) 



and, if is the hne obtained by contracting i^^j^-^ with ^(^/ j/), 



ip;i 



^n' = E E X ((j, ^). if, ^')) ^ t) . ^35) 



i=l i'=l 



The term / V*{dilj) e ^*(*) in Eq. flB[32p is the determinant of the (n — s + x (n — 
s + 1) matrix ^"^(t) (here 2n = "^l^i \Pi\), with elements G'jj,(t) = tj{f),j{f')Gfj', 
where fJ'^T,j{f)e{l,...,s} and tjj' = tn'{(j,j')) ■ ■ ■ tnHjj'))- 

J V*{d^) e-^*(*) = det ^^(t) . P36) 

Plugging this into Eq. (]B[32p finally gives Eq. (15 [3p . In order to complete the proof 
of the claims following Eq. (I5I3P we are left with proving the following Lemma. 

Lemma B.l dPT{t) is a normalized, positive and a-additive measure on the 
natural a-algehra of [0, 1]'^^-'^. Moreover there exists a set of unit vectors Uj G , 
j = 1, . . . , s, such that tjj/ = uj ■ uj/. 

Proof. Let us denote by bk the number of lines i & T exiting from the points 
0' J ^ -^k, such that i ~ Xk- Let us consider the integral 



/ dti . . . / dts_i Y\ {tn'{e) ■ ■ ■ tn{i)-l) — 1 
,...,Xs-i ° ° eeT 



P37) 



compatible with T 



and note that, by construction, the parameter tk inside the integral in the l.h.s. 
appears at the power bk — I- In fact any line intersecting dXk contributes by a 
factor tk, except for the line connecting Xk with the point in Xk+i\Xk. See Fig. 

El 

Then 

s-l 

II {tn'ii) . . . tn(i)-l) = n tk'' ' (038) 
e&T k=l 

and in Eq. (IB[37p the s — 1 integrations are independent. One has: 



/ dti... dt,_i Yi (tn'ii) . . . tnw-i) = n / ^^'^ ^k'' = n ?r ' (039) 

Jo Jo f^^^ \J0 J fc^i Ofc 
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which is well defined, since 6fc > 1. Moreover we can write: 

E = E E ■■■ E . »o) 

X2,---,Xs-l X2 Xs-1 

compatible with T Xi fixed Xi,X2 fixed Xi,...,Xs-2 fixed 

where the *'s over the sums remind that all the summations are subject to the 
constraint that the subsets Xi, X2, . . . , must be compatible with the structure 
of T. Now, the number of terms in the sum over X^, once that T and the sets 
Xi, . . . , Xfc_i are fixed, is exactly bk-i, so that 

1=61. ..6.-2, (041) 

X2,---,Xs-l 
compatible with T 

and, recalling that bs-i = 1, 

J2 dt,... dt,_i n (t„,(,) . . . t„(,)_i) = 11^^' (0^2) 

X2,-,Xs-i "^^ "^^ teT k=l ^ 

compatible with T 

yielding to / dPrit) = 1. The positivity and cr-addivity of dPrit) is obvious by 
definition. 
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We are left with proving that, for any given sequence of subsets Xi C X2 C 
■ ■ ■ G Xs compatible with T, we can find unit vectors uj G M'^ such that tjjr = 
Uj ■ Uji. With no losso of generahty, we can assume that Xi = {!}, X2 = {1,2}, 
. . ., Xs_i = {1, . . . ,s — 1}. We introduce a family of unit vectors in R'^ defined 
as follows: 



Vl 



-J = tj-i'^j-i + v^yi - , J = 2, . . . , s , 

where {vj}!^]^ is an orthonormal basis. From this equation, as desired: 



(E143) 



uj-uj, = tj...tj,_i. P44) 
This concludes the proof of the determinant formula. ■ 



Appendix C: Symmetry properties 

In this Appendix we prove Lemma El 2. The key remark underlying the nice 
properties stated in the Lemma is that both the Gaussian integration PM^dip) 
and the interaction V{ip) are invariant under the action of a number of remarkable 
symmetry transformations, which are also preserved by the multiscale integration. 
In the following, we denote by cxi, cr2, cts the standard Pauli matrices: 






Moreover, in order to avoid confusion with the Pauli matrices, we shall use the 
symbol r as spin index rather than a. As a preliminary result, let us start by 
listing all the transformation properties under which our theory is invariant. 

Lemma C.l For any choice of M,P,L, the fermionic Gaussian integration 
P{dijj) and the interaction V{ip) are separately invariant under the following 
transformations: 



(1) Spin flip: 'i/'k r ^ V'f 



k,— T' 



(2) Global U{1): ^-^^^ e*^"-?/'^^^, with e M/27rZ independent ofk; 

(3) Spin S0(2): | t^^^^' ] e^^^'^^ [ ^k,t,- | ^ ^^^j^ q eR/2nZ independent ofk; 

(4) Discrete rotations : i^^^ ^ e^'-'^^^'^^^'^^^'ip:^^^ and ip^^ — )• ■i/'^,^ ^e*'^^^^"''^)^, 
withTk= {ko,e~'^''4); 
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(5) Complex conjugation: V'kr V'-kr ^'^^ ^ ~^ ^* ! where c is generic constant 
appearing in P{d}^) or in V(\l/); 

(6. a) Horizontal reflections: ip^^ ^^'^~Rh^T '^^'^ "^kr ~^ "^^hkr^i' '^'^^^ ~ 
(6.h) Vertical reflections: ^/'^r V'|j„kT' — (^o, ^i, —^2); 



(^7j Particle-hole : V^kr ^"V^Pkr' "^kr ^V^Pkr' '"^^^^ ~ (^0; ~^l5 



2;; 



("5^ Inversion : ip^^^^ -ifTs'tpjk^r' V^k.r "'^'^/k.r'^S; /k = (-/cq, ^1, ^2)- 
Proof of Lemma IC.ll Let us first recall the definitions of PM{dip) and V{ip): 

PMm = j^ n #k,.,i^4,.,i*,.,2#k,.,2e-^<^^'^''''^-, 



where A/" is the normalization constant in Eq.f l4[8p and n± = (1 ± (T3)/2. Given 
these definitions and recalling the definition of ^j^^, we see that the invariance of 
PM{dip) and V{ip) are equivalent to the invariance of the following combinations: 

{**) ■= J2 fep.t'^^^k.t) (^k'-p,4'^"^k',J • (03) 

k,k',p 

Now, the invariance of (*) and (**) under symmetries (1) and (2) is completely 
apparent. Let us check one by one the invariance under the other symmetries. 

Symmetry (3). Note that (*) and (**) can be rewritten as 

(*) = E E v^k,^,p4~,.,p' , 04) 



k,p,p' 



I Yl ^^''^PP E fep,r,p4,r,p) (V^k'-p.r'.pV'k'.r'.p) 05) 



2 

k,k',p t,t' 



Invariance of these expression under symmetry (3) follows from the invariance of 
the combination Xlr ^k r p^k' rp'- fs^'Ct, denoting by the matrix e"*^'^^, we 
see that under (3) 
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which is invariant, simply by the fact that is orthogonal. 
Symmetry (4)- Under the action of (4), 

^ \voVL[k) tko J 

with Tk = (fco, R-fk). The r.h.s. of Eg.flClTl) can be rewritten as 

which is the same as (*), as it follows by the remark that i7(A;)e~*^^''^^'^^'' = 

27r ~* 

il(R~k). Regarding the interaction term, note that 

which immediately shows that (**) is invariant under (4). 
Symmetry (5). The term (*) is changed under (5) as: 



k 



which is the same as (*) (simply because VL*{k) = Q{—k)). Similarly, the term 
(**) is changed under (5) as: 



k,k',p 



which is the same as 



Symmetry (6. a). The term (*) is changed under (6. a) as: 



, . 7+ I iko VoVL*{k) \ 



k,T 



k,T 



Using the fact that VL*{k) = Q{{—ki,k2)), we see that this term is invariant 
under (6. a). Regarding the interaction term, if we note that aiUaCri = n_a, we 
immediately see that it is invariant under (6. a). 
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Symmetry (6.h). Invariance of the term (*) follows from the remark that 
VL*{k) = Q{{ki, — ^2)); invariance of (**) is trivial. 

Symmetry (7). The term (*) is changed under (7) as: 



5:*^-(„„:<;Vt?')'^*-' 013, 



which is the same as (*) (simply because Q*{k) = Q{—k)). Similarly, using the 
fact that = n^, the term (**) is changed under (7) as: 

a=± 

fe,t^°^P(k+p),t) fe',4^«^P(k'-p),J ' (014) 

k,k',p 

which is the same as (**). 

Symmetry (8). The term (*) is changed under (8) as: 

^ \von{k) tko J 



k,T 



which is the same as (*). Moreover, note that crsriocrs = ria, so that the term 
(**) is obviously invariant under (8). This concludes the proof of Lemma [C. II ■ 

Now, we are ready to prove Lemma [6l2. 

Proof of Lemma [6l2. The key remark is that, in addition to the fact that 
PM{dip) is invariant under the transformations (l)-(8), the ultraviolet and in- 
frared integrations P^dip^^'""'^) and P{d'tp^^'^''^) are separately invariant under the 
analogous transformations of the fields and ■ip^'-'^'\ Therefore, the effec- 

tive potential Vo('?/'^*'''''*) is also invariant under the same transformations. Let us 
restrict our attention to the quadratic contribution to Vo, i.e., 

P16) 

which, as we said, must be invariant under the symmetries (l)-(8) listed in Lemma 
IC.ll The very fact that we can write the quadratic term in the form Eq.( IC[16|l . 
with M^2(k) independent of the spin index r, is a consequence of symmetries (1)- 
(3). It is straightforward to check that Eq.f lGIll]) follows from the invariance of 
Eq. flC[l6l) under (4)-(8). 
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So, we are left with proving Eq.f l6ll2p . Let us start by showing that W{pp) = 
0. Writing W{p%) = a^I + a^ai + 020-2 + 030-3 and using the fact that, by 
Eq.f l6TTT|) . ly(p'J^) = aiW{p''p)ai = -a3W{p'^)cr3, we immediately see that 

= = = 0. Moreover, using the fact that (still by Eq.f l6TTT|) ) W{p'^) = 
e''^^^W{p'^)e-'^^^, we get a^cii = a^(cos(27r/3)o-i - w sin(27r/3)(T2), which 
implies = 0. 

Let us now look at A"^ := dk^Wip"^). Writing A'^ = b^I + + h%a2 + 63^3 
and using that, by Eq.f l6TTT|) . = -(Ag'^)* = A^'^ = aiA'^ai = a^A'^as, we 
see that Aq = —izol, with zq real and independent of u. In a similar way, 
writing A'^ = cp + d(ai + c^o-2 + c^o-3 and A% = dp + rf^o-i + rf^o-2 + rf^o-3, 
using the fact that = -(A^^)* = = -aiA^ai = -{A^'^)'^ = -a^A'^as 
and = -(^2"")* = -A^"^ = oiA%oi = -{A^'^f = -a^A'^a^, we see that 
Ai = C2U2 and A2 = udiai, with C2 and di real and independent of u. Finally, 
using the fact that, again by Eq- flHTTTD . = J2'^j^^e''^^''^{R^)ijA'^ , we get 
C2 = di =: —60 that, if combined with our previous findings, implies Eq.f l6[12p . 
This concludes the proof of Lemma [61 2. ■ 
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